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13.1  FlKDAMEaTTALS  OF  FEBffiACK  CXHTRDL  THECKY 


Pilots  might  be  inclined  to  associate  the  phrase  "control  system"  with 
only  aircraft  fli^t  control  systems.  Althouc^ti  the  control  system  theory  of 
this  course  has  a  large  ^plication  to  fli^t  control  systems,  this  material 
applies  to  any  process  or  ^stem  in  vAiich  control  is  exercised  over  some 
output  variable.  Examples  of  these  controlled  variables  are:  the  speed  of  an 
automobile,  the  tenperature  of  a  room,  the  attitude  of  a  spacecraft,  ad 
infinitum. 

Feedback  control  system  theory  is  often  called  several  different  things. 
It  mi^t  be  found  under  any  of  the  following  headings  or  titles:  Oontrol 

terns,  Servo-Mechanisms,  or  our  term.  Feedback 


Automatic  Control 


antroJ 


First,  the  difference  between  "open-loop"  and  "closed-loop"  control  will 
be  discussed.  Consider  the  roll  channel  of  an  aircraft  fli^t  control  system 
in  which  the  pilot  iiput  is  assumed  to  be  a  rate  ccrtmand.  That  is,  the  pilot 
ooranands  a  roll  rate  (J)  proportional  to  the  stick  displacement.  Figure 
13.1  shows  a  diagram  of  this  system.  Ihe  input  is  a  low  power  irput 
representing  a  selected  .value  of  roll  rate. 


POWER 
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ISlis  ittechanical  si9nal  is  tiien  anpllfied  in  a  I'Q'draulic  control  valve/actuator 
conbination  to  position  tli©  ailerons  accordingly  (5^)  •  deflected  ailerons 

then  react  with  the  airstream  to  produce  a  roll  nonetit  in  the  required 
direction.  The  nagnitude  of  the  resultant  roll  rate  Is  a  function 

primarily  of  the  dynamic  pressure  (q)  and  the  roaroent  of  inertia  about  the 
longitcdinal  axis  (I^) .  Both  are  part  of  aircraft  c^amics. 

There  are  inramerable  exanples  of  the  "non-feedback"  or  "open-loop"  type 
of  control  ^sterns.  For  instance,  a  gasoline  engine  in  an  automobile  has  a 
low  power  ii^nit,  the  throttle  position,  which  controls  the  speed  of  the 
vehicle  ejqjending  a  large  anoint  of  power.  In  a  simple  electronic  anplifier  a 
very  input  signal  controls  a  much  larger  output  signal.  In  all 

opei^loop  oontr^  systems  the  output  has  no  influence  on  the  iiput  vhatsoever. 
The  iiput  quantity  controls  the  output  only  directly  throu^  the  intermediate 
oonponents.  Referring  back  to  the  eiircraft  roll  control  system,  a  lateral 
stick  displacement  of  a  specified  amount  will  not  command  a  constant  roll  rate 
all  conditions.  As  the  conditions  in  the  intermediate  coiponents 
change,  such  ais  the  dynamic  pressure,  moment  of  inertia,  hydraulic  pressure, 
taiperature  of  l^^aulic  fluid,  condition  of  hydraulic  coiponents,  temperature 
effects  on  modulus  of  elasticity  of  metal  components,  etc.,  the  resultant  roll 
rate  for  a  specified  iiput  will  vary. 

The  performance  of  any  control  ^stesm  with  respect  to  maintaining  the 
output  quantity  as  close  as  possible  to  the  input  quantity  can  be 
substantially  improved  by  feeding  back  the  output  for  ccnparison  with  the 
iiput.  The  use  of  the  difference  resulting  from  this  ccnparison  as  an 
actuating  signal  constitutes  a  feedback  or  closed-loop  control  system. 


FIGURE  13.2.  CDOaED-KXDP  CCNTRDL  SYSTEM 
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Figijre  13.2  shows  how  the  open-loop  system  of  Figure  13.1  can  be  changed  to  a 
closed-loop  ^tera  by  the  addition  of  an  outer  feedback  loop  to  corpare  the 
input  vd.th  the  output.  Thus,  the  effect  of  variations  in  the  intermediate 
conponents  can  be  eliminated  in  that  a  corrective  signal  (e)  will  continue  to 
exist  xmtil  the  output  properly  matches  the  input. 

A  serious  disadvantage  of  closed-loop  control  systems,  however,  is  that 
they  can  make  an  otherwise  stable  system  unstable.  Ihe  possibilily  of 
instability  is  the  prime  reason  for  the  existence  of  the  science  of  feedback 
control  systan  analysis.  Uie  first  and  major  effort  in  control  system 
analysis  is  the  determination  of  vdiether  or  not  the  closed-loop  systan  is 
stable.  After  this  fact  is  established,  other  response  characteristics  may  be 
found. 

Stability,  with  respect  to  control  systems  is  defined  as  follcws:  A 
stable  system  is  a  systan  in  vAiich  the  transients  die  out  with  increasing 
time. 

13.2  NQMEI^CIATUPE 


The  following  nonenclatxrre  is  used  in  this  chapter: 

R  =  input  variable 
C  =  output  variable 

Each  of  these  might  represent  aiy  quantity  depending  on  the  system  such  as 
angular  or  linear  position,  current,  voltage,  degrees  of  temperature,  etc,  or 
the  time  rate  of  change  of  those  above. 

Ihe  following  symbol  represents  a  suirmer  or  differential.  It  indicates 
the  algebraic  sumnfiation  of  the  input  quantities  according  to  the  arrows  and 


t- 


FIGDRE  13.3.  SIMMER  OR  DIFFERENTIRL 
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the  signs.  Ihe  exaiiple  in  Figure  13,3  shews 

e  =  R-C  (13.1) 

The  symbol  for  gain  or  aitplification  factor  is  K. 

Control  ^tems  are  generally  described  through  the  \jse  of  block  diagrams 
as  in  Figures  13.1  and  13.2.  However,  instead  of  words  to  indicate  the 
process  or  aeration  occurring  witiiin  each  block,  there  appears  v^at  is  called 
a  transfer  funccion  (Figure  13.4).  The  term  "transfer  function"  might  be 
thou^t  of  as  vAiat  is  done  to  the  input  to  produce  the  output.  Althou^  the 
transfer  functions  within  the  blocks  are  generally  written  in  terms  of  seme 
operator  notation,  th^  are  often  described  graphically,  especially  for 
nonlinear  systems.  A  definition  of  transfer  function  is:  the  ratio  of  the 
output  to  the  input  expressed  in  Laplace  operator  notation,  assuming  zero 
initial  conditions.  Ihe  transfer  function  is  essentially  a  mathematical  model 
of  the  systan  and  entoodies  all  the  physical  characteristics  of  the  system 
i.e. ,  mass,  dairping,  etc. 
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FIGURE  13.4.  TRANSFER  FUNCTION 
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13.3  DIFFERENTIAL  BQUATIOt^  -  CLftSSlCAL  SOLUTIONS 

Differenti2a  equations  for  a  control  system  will  illustrate  the  types  of 
responses  to  be  expected  from  first  and  second-order  systems.  Biese  two 
exanples  are  used  throughout  this  course  because  higher-order  systems  produce 
a  transient  response  consisting  of  the  sum  of  first  and  second— order 
responses. 

The  reason  the  transient  response  is  significant  rather  than  the  steady 
state  or  complete  response,  concerns  the  stability  of  the  system.  Since 
positive  stability  requires  that  the  transients  die  away  with  increasing  time, 
the  transient  solution  of  the  differential  equation  describing  the  system  is 
most  i]({X>rtant  to  the  analysis.  The  transient  solution  also  provides  other 
inportant  response  characteristics. 

Figure  13.5  shows  a  sinplified  block  diagram  of  VTOL  Auto  Pitch  Control 
with  inertia,  1^,  and  limited  aero-danping  proportional  to  pitch  rate.  Pitdi 
attitude  is  maintained  by  reaction  control  jets.  These  jets  produce  a  torque 
proportional  to  a  valve  position.  Torque  -  |/e,  where  /#  is  the  gain  of  the 
valve  and  e  the  input  to  the  valve.  The  loop  is  closed  ly  conparing  the 
output  pitch  attitude  to  the  commanded  pitch  attitude. 


COMMANDED  OUTPUT 


FIGURE  13.5  VTOL  AUTOMATIC  PITCH  CCmROL  BLOCK  DIACSEy^M 


OMs  cotparison  produces  the  error  signal,  e,  viiich  is  the  iiput  to  the  linear 
valve.  Bie  resulting  torque,  pe,  is  applied  to  the  v^iicle  to  diange  the 

pitch  attitude. 

TWO  situatiOTS  will  be  considered  in  order  to  sinplify  the  problem.  In 
the  first  case,  only  the  effect  of  viscous  danping  will  be  considered.  This 
will  result  in  a  first-order  differential  equation.  The  second  case  will 
both  inertia  and  viscous  friction  and  will  result  in  a  second-order 
system.  The  first  aid  second-order  differential  equations  will  be  solved  for 
the  transient  response. 

13.3.1  First-Order  Systan 

Using  Figure  13.5,  including  only  the  effect  of  danping  on  the  vehicle, 
the  differential  equation  of  the  system  can  be  written  by  equating  the  ^plied 
torque  to  the  absorbed  torque.  The  torque  applied  by  the  reaction  jets  is 
by  the  viscous  friction  (aero-danping)  of  the  v^cle. 

pe  =  be^ 

The  output  of  the  conparator,  e  =  ”  ®0  ^stem  differ¬ 

ential  equation  or  the  equation  of  motion 

p(e^-eo)  =  be^ 

ipplied  Torque  =  Absorbed  Torque 

Using  the  operator  "p"  notaticm  (vdiere  p  =  d/dt)  to  determine  the  system 
transient  response,  the  homogeneous  equation  becomes 

7  P  ®o  ®o  = 


0 


'A 


the  root  of  the  characteristic  equation  is 


P  » 


Since  the  transient  response  is  assumed  to  be 


00  (t) 


transient 


n  p  t 
I  C,e 
i-1 


(13.2) 


vdiere  p^  »s  are  the  roots  of  the  characteristic  equation,  the 
t^®^sient  response  for  our  first— order  system  is 


©oft) 


transient 


(13.3) 


For  positive  gain  //  and  damping  factor  b,  ej(t)is  always  stable.  Bius,  a 
first^rder  system  has  only  real  roots  of  the  characteristic  equation  and  the 
transient  response  is  either  an  exponential  increase  or  decrease  depending  on 
the  sign  on  the  time  constant.  The  time  constant  is  the  reciprocal  of  the 
coefficient  of  t  in  the  exponent  of  e  (b/t/  in  our  case). 

The  time  constant,  generally  given  the  symbol  t,  can  be  defined  as  the 
value  of  time  that  makes  the  exponent  of  e  equal  to  -1.  m  one  time  constant 
the  exponential  has  decreased  from  the  value  1  to  the  value  .368.  Figure 
13.6  shows  a  plot  of  the  transient  response  of  a  first-order  stable  system. 
Time  constants  are  discussed  in  more  detail  in  Paragraph  13.5.4. 
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FEURE  13.6.  PliDT  OF  EERST  ORDER  TRSNSIENT  RESPCXJSE 


A  stable  transient  response  requires  the  root  of  the  characteristic  equation 
be  negative. 


13.3.2  Second-Order 


Equating  applied  torque  to  absoi±)ed  torque  of  Figure  13.5  again,  but 
including  the  inertia  and  daitping  effects  yields 


ve  =  be.  +  le. 


again 

6  =  e.  -e^ 

and 

^(©1  -  Sq)  =  bOp  +  I  'e'^ 
iOq  +  be^  +  ye^j  =  pe^ 

is  the  equation  of  notion  of  the  system.  Bie  hcmogeneous  equation  in  operator 
notation  is 
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0 


+  +  y)  = 


Hie  characteristic  eguation  is 

Ip^  +  bp  +  y  =  0 


vAiose  roots  are 


-b  ±Vb^  - 


■1.2 


4y  I 


21 


(13.4) 


(13.5) 


Dspending  on  the  relative  magnitudes  of  the  gain,  dairping  factor,  and 
inertia,  the  roots  of  the  characteristic  equation  mi^t  be  real  or  ccriplex 
ther^y  indicating  different  types  of  response.  If  the  roots  turn  out  to  be 
real,  the  transient  re^onse  is  merely  the  sum  of  the  two  resulting  first- 
order  ejqponential  terms.  If  the  roots  are  catplex,  however,  they  always 
appear  in  ccnplex  conjugate  pairs  in  the  following  form; 

Pl,2  =  o  +  jw^ 

vhere  o  is  Ihe  real  part  and  the  imaginary  part  of  the  roots.  Hiese 
ccnplex  roots  yield  a  solution  of  the  form 


®  0  ^  transient  ^1® 


(a+jto^)t  (a  -  jto^)t 


+ 


(13.6) 


After  ccnplex  variable  manipulations,  this  eiqpression  can  be  shown  to  be 
equivalent  to 


eQ(t)  =  Ae'^^  cos((i)^t  +  (}>) 


(13.7) 


vhere  A  and  (|)  are  derived  from  the  coefficients  and  C2. 
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This  is  the  form  of  the  solution  viienever  the  characteristic  equation  has 
coiplex  conjugate  roots. 


FIGURE  13.7.  EXPONENTIALLY  DAMPED  SINUSOID  - 

TYPIC3VL  SECXM)-C®DER  SYSTEM  RESPONSE 

It  is  called  an  ejqponentdally  damped  sinusoid  and  consists  of  a  sine  wave  of 
frequency  oi^  vihose  magnitude  is  Ae*^;  that  is,  it  is  decreasing  esqponentially 
with  time  if  o  is  a  negative  quantity.  A  typical  second-order  response  is 
plotted  in  ligure  13.7. 

Inferring  back  to  the  solution  of  the  characteristic  equation.  Equation 
13.5,  the  real  part  can  be  recognized  as  the  exponent  of  e 
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k 


(43.8) 


frequency  of  the  oscillation  of  the  transient 


».  -  (13.9) 

V  41^ 

Uie  quantity  b  represents  the  effective  damping  of  the  system,  if  b  equals 

2  TJJT  the  two  roots  p^  ^  are  equal.  This  is  the  critical  level  of  damping  and 
is  written  b'  -  2 

Bie  damping  ratio  is  defined  as  the  ratio  of  actual  damping  to  the 
critical  valiie  of  danping 


r  .  actual  cteyinq  b  b 

critic2d  danping  *  * 

2 


(13.10) 


When  C  IS  greater  than  zero  but  less  than  one  the  roots  are  complex  and 
the  solution  is  a  danped  sinusoid  of  the  form  of  Figure  13.8  and  is  called 
tmderdamped.  When  C  is  greater  than  one  the  roots  are  real  and  the  response  is 
overdaoped.  Vlhen  C  is  negative,  the  system  is  unstable. 

Bie  undanped  natural  frequency,  is  defined  as  the  frequency  of 
oscillation  of  the  transient  if  the  danping  is  zero.  Prom  Equation  13.9 


\/^'  (13.11) 


Bie  response  in  the  case  of  no  danping  is  a  sine  wave  of  constant  anplitude. 

Second-order  equations  (or  factors  in  more  complex  systems)  are 
frequently  written  in  terms  of  the  danping  ratio  and  undamped  natural 
frequency.  Factoring  ^  from  Equation  13.4  leaves 


+  ^P  +  1  “  0  (13.12) 


13.11 


vrtiere  I///  can  be  recognized  as  and  h/ti  equals  2C/\‘ 

equation  becomes 


The  characteristic 


JL  p=  +  25  p  +  1 

n  n 


0 


Multiplying  by  produces  the  standard  form  of  the  second-order  system. 


+  2Z;»„p  +  =  0 


(13.13) 


The  roots  of  this  equation  are 


p  =  a  ± 


-  ±  j«>„  >/l  -  ^ 


(13.14) 


And  the  transient  response  in  terms  of  C  and  is 


c(t) 


transient 


.  — Cco  t 
Ae 


cos 


(13.15) 


Figure  13.8  shows  a  family  of  curves  representing  the  response  to  a  step  input 
of  a  second-order  system  as  a  function  of  C.  These  curves  illustrate  the  fact 
that  the  amount  of  overshoot  and  the  time  to  arrive  at  the  input  value  are  a 
function  of  C« 
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r' 


FIGURE  13.8.  SECQLnD-ORDER  TRANSIEtCT  RESPC»ISE  VERSUS  ? 
13.4  TRANSFER  FUNCTIONS 


Ihe  first  and  second-order  equations  derived  in  the  previous  paragraph 
were  solved  tising  the  classical  method  to  show  the  types  of  response  to  be 
expected  fran  each  type  system.  In  practice,  this  approach  is  extremely 
laborious,  if  not  impossible,  for  more  ccmplicated  systems.  Bierefore,  more 
advanced  techniques  are  used  vhich  do  not  produce  the  total  solution  but  do 
indicate  \(hether  or  not  the  system  is  stable;  and  if  not,  provide  inficrmation 
about  how  to  make  the  system  stable.  These  sophisticated  techniques  generally 
use  Laplace  transforms.  The  use  of  operational  calculus  offers  a  definite 
advantage  in  that  transfer  functions  can  be  manipulated  using  the  normal  rules 
of  algebra.  It  also  iitposes  a  severe  restriction.  The  systems  to  be 
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analyzed  must  be  representable  by  linear  differratial  equations  with  constant 
coefficients.  The  method  of  determining  these  transfer  functions  will  now  be 

described. 

First,  we  will  consider  the  ^tem  of  Figure  13.5  in  vdiich  only  danping 
was  included.  Its  equation  of  motion  is 


letting  t  =  b/p  we  have 


tGo  +00=6^ 


TaJdng  the  Laplace  transform  using  the  notation  ^[e^Ct)]  -  6q(s) 

ts0q(s)  -  xegCO^)  +  Bq  (s)  =  Gj^Cs) 


where 


0^(0"*^)  is  the  value  of  6Q(t)  at  t  =  0^ 

00  (s)  [ts  +  1]  =  e^(s)  +  tGq  (o’*") 


0^(s)  tGq  (0^) 

~  TS  +  1  TS  +  1 


Thus  we  see  that  the  input  to  the  systan  is  acted  ipon  by  the  transfer 
function 

^  (13.16) 

TS  +  1 

and  also  the  initial  condition  is  acted  upon  by  this  transfer  function. 
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Fran  this  brief  discussion  we  can  see  that  if  we  assume  all  the  initial 
conditions  to  be  zero  we  obtain  the  relationship 


00  (s) 


1 

TS  +  1 


(13.17) 


Our  first-order  system  can  then  be  described  in  the  manner  of  Figure  13.4 
vAiere  G(s)  =  1/(ts  +  1),  the  transfer  function. 

The  transfer  function  of  our  second-order  system  in  vdiich  the  inertia  and 
daitping  vere  considered  will  now  be  determined  using  the  same  procedure.  From 
its  equation  of  motion, 


1  Oj,  +  b  e„  +  peg  =  pe^ 

taking  the  Laplace  transform  and  assuming  all  initial  conditions  to  be  zero  we 
have 


eo(s)  eg(s)  +  6o(s)  =  e^(s) 

The  transfer  femction  is  then 


00  (s) 

0^(3) 


is2 

y 


+  1 


(13.18) 


The  transfer  functions  that  have  been  developed  for  first  and  second- 
order  systems  (Equations  13.17  and  13.18)  are  obtained  frcm  the  equation  of 
motion  of  the  v^ole  system  with  the  feedback  loops  closed  (Figure  13.5). 
Therefore,  th^  are  called  closed-loop  transfer  functions. 

The  denominator  of  Equation  13.18  is  equivalent  to  Equation  13.12  and  is 
the  characteristic  equation  of  the  system. 
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13.5  TIME  DOM?mi  ANALYSIS 


Much  of  the  work  of  the  control  system  engineer  is  done  in  the  s-donain 
to  take  advantage  of  simplicity  of  solution,  but  the  response  of  a  system  is 
in  the  domain.  Die  time  response  of  a  ^stem  is  divided  into  two  parts: 
(1)  the  transient  response,  and  (2)  the  steacfy-state  re^nse. 

<=<«  -  'trans^*  's.s.  <« 


In  order  to  analyze  a  control  system,  we  discuss  the  performance  of  the  systan 
in  terms  of  time  response  to  a  specific  input.  Bor  a  given  ^stem,  a  specific 
ii^jut  will  result  in  a  predictable  transient  response  and  a  steacty-state 
error.  Control  qrstem  performance  specification  can  be  stated  in  terms  of  the 
transient  bshavior  of  the  system  and  the  allowable  stea^— state  error.  In 
general,  the  steai^-state  error  can  be  a  fimction  of  time;  however,  we  usually 
vrant  lim  e(t)  =  0. 

In  reality,  control  system  specification  and  obtainable  real  world 
solutions  are  a  oonpromise.  Die  first  order  of  business  in  analyzing  a 
control  system  is  to  determine  if  the  s^^stem  is  stable.  If  it  is  stable,  then 
it  will  be  tested  to  determine  if  it  meets  the  performance  specifications. 
Die  response  of  the  system  to  specific  test  inputs  will  provide  several 
measurements  of  performance. 


13.5.1  Typical  Time  Domain  Test  Input  Signals 

13.5.1.1  Die  step  input  is  the  most  ccauronly  used  test  signal.  Diis  input  is 
simply  an  instantaneous  change  in  the  reference  input  variable  (Figure  13.9) . 


r(t) 


R 


u_i(t) 


O 


t 


l(t)  =  R,  t  >  O 
r(t)  =  0,t<0 
R  =  CONSTANT 


FIGURE  13.9.  STEP  INPUT 
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r(t)  =  u_j(t) 

vrtiere  u_^(t)  is  tte  tmit  step  function.  The  quantity.  r(t)  is  not  defined  at 
t  =  0.  The  laplace  transform  of  the  unit  step  is 

=  I  (13.20) 


Therefore  the  Laplace  of  r(t)  =  Ru_j^(t)  is  R/s. 

13.5.1.2  Ranp  Function.  The  rartp  signal  is  the  integral  of  the  unit  step  and 
is  often  called  the  velocity  input  (Figure  13.10) . 


FIGURE  13.10.  RAI’IP  INPUT 


r(t)  =  Rt  u_^(t) 


and 


oZfCRt  u^^(t)3  =  ^ 


(13.21) 


13.5.1.3  Parabolic  Input.  The  parabolic  input  signal  (Figure  13.11)  is  the 
integral  of  the  rartp  signal  and  is  often  referred  to  as  the  acceleration 
input. 
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FIGURE  13.11.  PARABOLIC  INPUT 
r(t)  =  Rt^u_^(t) 


and 


(13.22) 


13.5.1.4  Power  Series  Input.  An  ii?jut  made  up  of  the  sums  of  a  step,  ranp 
and  a  parabola  would  be  a  power  series  of  pcwer  2. 


r(t)  =  R 


r  t^l 

l  +  t+|-  u_^ 

ai 


(t) 


(13.23) 


13.5.1.5  Unit  Iitpulse.  Another  useful  iiput  is  the  unit  impulse  (Figure 
13.12). 


»0,  ELSEWHERE 


FIGURE  13.12.  UNIT  IMPIESE 
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As  e  -*■  0,  the  function  r(t)  approaches  the  inptilse  function  6  (t) . 

c^C6(t)3  =  1  (13.24) 

13.5.2  Time  Response  of  a  Second-Order  System 

Consider  the  closed  loop  block  diagram  in  Figure  13.13 


-  ♦ 


^  C(s) 


FIGURE  13.13.  CLOSED-LOOP  CONTROL  SYSTEM 


vAiere 


C(s) 


^  G(s)  R(s) 
1  +  G(s) 


Let 


G(s)  — 


K 


s(s  +  a) 


K 


C(s)  = 


s(s  +  a) 


R(s) 


1  + 


K 


s(s  +  a) 


R{s) 


s  +  as  +  K 


This  equation  can  be  generalized  in  terms  of  c  and  u 


n 


(13.25) 
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Let  a 


2C(0„r  K  *  * 


Ihen  the  control  ratio  C(s)/R(s)  is 


C(s) 

R(s) 


+  20a.  s  +  <0.  ^ 

n  n 


(13.26) 


For  a  unit  step  ii^t,  R(s)  =  1/s,  and 


“  iHr 


s  (s  +  20*).  s  +  to.  ) 


(13.27) 


Taking  the  inverse  Laplace  gives  the  transient  response 


c(t)  »  1-  ,  e  “^"n^sin  («  tVl-  +  ♦) 


(13.28) 


v^ere  ♦  -  tan"^  Vl  -  ^ The  transient  response  of  this 

system  varies  according  to  the  selected  value  of  C»  Figure  13.14  depicts  this 

variation. 

Several  standard  performance  specification  terms  common  throughout 
industry  are  illxistrated  in  Figure  13.15. 
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FIGURE  13.15.  TIME  DOMAIN  SPBCIFICAnC»IS 


Biese  transient  performance  specifications  are  usually  defined  for  a  unit 
step  ii^t. 

Oyershoot  -  indicated  by  largest  error  between  input  and 
output  &ring  the  transient  state.  We  can  determine  the 
magnitude  by  using  the  previously  developed  equation  for 
a  step  input  to  a  second— order  system. 

-  1  _  ^  e  sin  ^w^t  \/  1  -  +  tan~^  ~  7  ^  )(13.29) 

^1  -  ^  ^ 

Taking  the  derivative  of  this  equation  and  equating  to  zero  yields 

-  0  -  —  e"^“ntsin  /  «  t\/l  -  C"  +  tan'^  ~  ^  \ 

/l  -  ^  / 


1 _ 

\/l  -  ^ 

This  derivative  is  zero  \dien  <o^t  \/l  -  =  0,  n,  2n,  etc. 

The  peak  overshoot  occurs  at  the  first  value  after  zero  (with  initial 
conditions  equal  to  zero). 

Therefore,  the  time  to  maximum  or  peak  overshoot  is 


1. 


c{t) 


(13.31) 


Substituting  this  value  into  c(t)  yields  the  peak  response,  M  . 


M 


".v  1  -  C' 
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Hie  overshoot  for  the  unit  step  ii^t  is 


Overshoot  ■  -  1  -  exp 

and  the  percent  of  overshoot 


Cn 


\/TT 


(13.33) 


P.O. 


”p  "  ^  X  100% 


100%  e3^ 


Cn 


v/TT 


(13.34) 


2. 


Delay  Timey  is  the  time  required  for  the  response  to  a  unit  step 
to  reach  50%  of  its  final  value. 


3.  Rise  Time^  is  the  time  required  for  the  response  to  a  unit  step 
to  rise  from  10%  to  90%  of  its  fined,  value. 


4.  Settling  Time,  Tg,  is  the  time  required  for  the  response  to  a  unit 
step  to  decrease  to  and  to  stay  vdthin  a  specific  percentage  of  its 
final  value.  Ctommonly  used  values  are  2%  or  5%  of  the  final  value. 

13.5.3  Hi^Ter-Order  Systems 

Ihe  relationships  developed  in  the  preceding  paragraph  using  (o^  and  c 

apply  equally  well  for  each  ccnplex-conjugate  pair  of  poles  of  an  n^-order 
system.  Ihe  distinction  is  that  the  doninant  ^  and  ajply  for  that  pair  of 
complex-conjugate  poles  vdiich  lie  closest  to  the  imaginary  axis.  The  values 
of  c  and  0)^  are  dominant  because  the  corresponding  transient  terra  has  the 
longest  settling  time  and  the  largest  magnitude.  Therefore,  the  dominant 
poles  primarily  determine  the  sh^je  of  the  time  response,  c(t).  A 
nondominant  pole(s)  has  a  real  axis  component  that  is  at  least  six  times 
further  to  the  left  than  the  corresponding  corponent  of  the  dominant  pole(s). 
Components  of  c(t)  due  to  nondcminant  pole(s)  die  out  relatively  quickly,  and 
can  be  neglected.  (13.1:245), 

13.5.4  Time  Constant,  t 

The  time  constant  is  used  as  a  measure  of  the  exponential  decay  of  a 
response.  Etor  first-order  systems,  the  transient  response  is  an  exponential 
function  described  by  Ae"™^,  Figure  13.16. 


FTEUPE  13.16.  PLOT  CF  EXPCNENTJAL  e"”^^ 
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Bie  value  of  time  v^ch  makes  the  ejqjonent  of  e  equal  to  -1  is  defined. as 
the  time  constant,  t 


■^T  =  -1 


T 


In  one  time  constant  the  exponential  e  will  decrease  frcin  a  value  of 
1.0  to  a  value  .368.  Table  13.1  shows  values  for  other  time  constants. 


TABLE  13.1 
TIME  CONSTANT  TABLE 


t 

e 

It 

.368 

2t 

.135 

3t 

.0498 

4t 

.0183 

Ohe  time  constant  is  another  way  of  specifying  settling  time.  The 
exponential  will  decay  to  37%  of  its  initial  value  in  t  seconds  (one  time 
constant) .  In  3t  the  exponential  is  within  approxdmately  5%  and  in  4t 
approxdmately  2%  of  its  final  value.  Rxr  a  second-order  imderdaitped  systan  of 
the  form 


c(t)  = 


sin  (cij^Bt  +  (|)) 


“infc 

the  response  is  bounded  by  the  exponent  of  the  form  (l/8)e  .  The 
Reifications  of  ?  and  determine  the  bounding  exponential  curve.  The  time 
constant,  t,  for  these  systems  is 


(13.35) 
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13.6  STABIIJIY  DETEIMINATION 


Q3i6  mosli  important  aroa  of  ths  anai^sis  of  a  closed~loop  control  systsni 
is  the  determination  of  stabilily.  A  system  is  said  to  be  stable  if  the 
out^t  of  the  system  corresponds  to  the  input  after  transients  die  out.  A 
system  is  said  to  be  unstable  if  the  transients  do  not  die  out  or  if  they  grow 
larger  following  a  disturbance. 

Stability  is  an  inherent  characteristic  of  the  :^stem  and  only 

i^n  the  system  itself,  not  tpon  the  input  or  forcing  function.  Hance,  if  a 
^tem  is  unstable,  any  input  will  cause  the  ^tem  to  diverge.  If  the  ^stem 
is  stable,  any  bomded  iiput  will  cause  a  boxmded  response. 

Bie  problem  in  determining  stability  is  ascertaining  vdiether  or  not  the 
transients  of  a  system  will  die  out  BEFORE  the  system  is  built. 

Ws  must  determine  the  conditions  under  vdiich  a  system  will 
unstable  and  be  able  to  tell  when  this  happens  in  the  analysis  of  the  system. 
Several  methods  are  available  for  determining  stability:  root  locus.  Bode 
plot,  Routh's  stability  criterion,  and  Nyquist  Criterion.  Only  the  root  locus 
and  the  Bode  plot  methods  will  be  presented  in  this  chapter. 

13.6.1  Stability  in  the  s-Plane 

Since  this  course  is  concerned  with  linear  systems,  i.e. ,  those  vhose 
differential  equations  i  are  linear  vidth  constant  coefficients,  the  transient 
response  is  of  the  form 


eft) 

'  'transient 


n  •  s.t 
=  Z  k.e^ 
i  =  1  ^ 


where  n  is  the  order  of  the  differential  equation  and  the  values  of  s  are  the 
roots  of  the  characteristic  equation  which  are,  in  general,  conplex. 


s  =  a  +  jo)^  (13.36) 

a  is  the  real  part  of  the  ccitplex  variable  s  and  is  the  imaginary  part  of 
the  ccmplex  variable  s.  Ihe  notation  used  to  indicate  this  is 
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0  =  Re  [s3 
(1)^  =  Hn  fs3 

Previously,  we  discussed  only  a  first  and  second-order  system  and  saw  the  type 
of  transient  response  to  be  es^jected  from  each.  Ihe  characteristic  equation 
of  hi^er-order  ^sterns,  however,  can  in  theory  be  factored  into  the  product 
of  several  first  and  second-order  factors  depending  on  the  order  of  the 
equation.  This  is  demonstrated  in  Equations  13.37  and  13.38 

+  .  .  .  .  Aq  =  0  (13.37) 

Can  be  ejqjressed 


0  (13.38) 


Bie  transient  re^onse  of  a  ocmplex  system  is  the  sum  of  those  associated  with 
each  of  the  first  and  second-order  factors.  Each  root  of  the  characteristic 
equation  must  be  of  one  of  the  forms  shown  in  Figure  13.17.  Opposite  the 
possible  values  of  the  roots  on  the  left  are  shown  the  corresponding  transient 
response  components  as  a  function  of  time.  Note  that  cctnplex  or  imaginary 
roots  always  occur  in  catplex-conj\:igate  pairs,  lhat  is,  they  have  imaginary 
parts  of  equal  magnitude  but  cure  opposite  in  sign. 

All  the  possible  values  of  s  can  also  be  described  throu^  use  of  a 
complex  plane  —  in  this  case  the  s-plane.  A  conplex  plane  is  one  in  vdiich 
the  value  of  the  real  part  of  the  conplex  variable  is  the  distance  along  the 
abscissa  and  the  magnitude  of  the  imaginary  part  is  described  along  the 
ordinate.  Ihese  are  called  the  real  and  imaginary  axes,  respectively.  'Qie 
crarplex  variable,  s,  is  then  a  position  vector  in  the  complex  s-plane  viiere 
o  =  Re  {s3  is  the  roagnitiade  of  the  real  cociponent  and  Cs3  is  the 


imaginary  component. 
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Figure  13.18  shows  the  s-plane.  If  the  values  of  s,  which  are  the  roots. of 
the  characteristic  eguation^  are  plotted  in  .the.  s~planef  they  produce  a 
transient  solution  cooponent  as  indicated.  Areas  in  the  s-plane  in  ^ich 
roots  produce  stable  and  unstable  responses  are  also  shown.  Roots  yielding  an 
undamped  response  or  neutrally  stable  output  are  all  on  the  imaginary  axis. 
Foots  associated  with  non-oscillatory  response  are  all  on  the  real  axis.  A 
root  of  the  characteristic  equation  at  the  origin  (s  ■=  0)  has  a  transient 
solution  equal  to  a  constant. 

The  mathematical  definition  of  a  stable  system  is  one  in  vdiich  the  roots 
of  the  equation  have  only  negative  real  parts.  In  other  words,  vdiere  s  = 

o  +  jiOj  are  the  roots  of  the  characteristic  equation,  <r  <  0  produces  a  stable 
system. 
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VALUE  OFs 


TIME  RESPONSE 


i 

PICTORIAL  SOLUTION 


a(POSITIVE) 


(7  (NEGATIVE) 


Ce®* 


±J“n 

( 4  zero) 


Ccos(«’„t  +  (i>) 


(a  POSITIVE) 


cos  { +  <t> ) 


FIGURE  13.17,  TRANSIENT  SOLUTION  CF  LINEAR 
CONSTANT  COEFFICIENT  BQUATICNS 
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C  cos  C«Jn  t  +  t) 


Ce»’t  COS  (’u«ft+^ 
(<r  NEGATIVE) 


^Ce«-‘  COS  (i^t+D) 
(<r  POSITIVE) 


STABLE 


UNSTABIB 


FIGURE  13.18  COMPLEX  s-PLANE 


3.3 
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13.6.2  Additional  Poles  emd  Zeros 

The  results  of  adding  a  real  pole  or  a  real  zero  to  the  basic 
second-order  control  ratio  as  given  ly  Equation  13.26  will  be  investigated. 
When  a  pair  of  complex— conjugate  poles  are  dominant/  the  approximations 
developed  in  Paragrajdi  13.5.2  yield  accurate  resxilts.  The  addition  of  a  third 
real  pole  to  a  seccaid-order  transfer  function  can  significantly  alter  the 
system  time  response  c(t)/  and  the  approximations  given  in  Paragraph  13.5.2 
no  Icmger  give  accurate  results. 

The  effects  of  a  third  real  pole  can  be  seen  by  considering  the  control 
ratio 


The  time  response  resulting  from  a  unit  step  input  is 


c(t)  ■=  1  +  2  |a(  e“^n^  sin  («^t  +  ♦)  +  Be^s*^  (13.40) 


(13.40b) 


The  transient  term  due  to  the  real  pole,  Pj ,  has  the  form  BePa*^,  where  B  is 
always  negative.  Therefore,  the  peak  overshoot  is  reduced,  the  settling  time, 
tg ,  may  be  increased  or  decreased,  and  the  phase  angle  is  decreased.  This  is 
the  typical  effect  of  adding  a  third  real  pole.  The  further  to  the  left  Pj 
is,  the  smaller  the  magnitude  of  B  and  the  effect  on  c(t).  Typical  i-nmo 
responses  as  a  function  of  the  real  pole  location  are  shown  in  Figure  13.19 
(13.1:350). 
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FIGURE  13.19  TIME  RESPONSE  AS  A  FUNCTIOJ 
OF  REAL  POLE  LOCATIC*^ 


The  time  response  is  also  altered  by  the  addition  of  a  real  zero  to 
second-order  transfer  function  of  Equation  13.26.  The  control 
becomes 


C(s) 

R(S) 


K(s 


(f  +  s  +  0)^^ 
The  time  response  resulting  from  a  unit  step  input  is 


1/2 


c(t)  e  1  +  ^  |(2  -  e  "  sin  (w^t  +  4.) 


\diere 


“d 


z  -  ttil 


the  basic 
ratio  now 

(13.41) 


(13.42) 


4»  «  tan' 


(13.43) 


Fran  Equations  13.42  and  13.43,  it  is  seen  that  the  addition  of  a  real  zero 
affects  both  the  magnitude  and  phase  of  the  transient  part  of  c (t) .  The  real 
zero  tends  to  increase  the  overshoot  and  decrease  the  phase  angle  of  c  (t) . 
This  effect  becones  more  dramatic  as  the  zero  approaches  the  imaginary  axis. 
This  is  illustrated  in  Figure  13.20  for  a  stable  second-order  system  with 
held  constant. 
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FIGURE  13.20.  TIME  RESPONSE  OF  A  SBCXM>-ORDER  SYSTEfl 
AS  A  FUNCTICN  OF  REAL  ZERO  LOCATION 

Several  things  can  be  seen  by  examining  the  responses  in  Figure  13.20.  First,- 
it  should  be  noted  that  rise  tire  is  decreased  and  the  overshoot  is  increased 
with  the  addition  of  a  real  zero.  Ihis  effect  is  more  noticeable  as  the  zero 
moves  closer  to  the  imaginary  axis.  This  is  to  be  expected,  because  when  the 
zero  is  at  the  origin,  it  acts  as  a  pure  differentiator  of  the  input. 
Differentiation  of  the  unit  step  input  yields  the  unit  irtpulse.  When  the  zero 
is  in  the  ri^t  half  plane,  the  response  is  stable  but  the  direction  of  the 


13.34 


initial  response  is  opposite  to  the  final  steady-state  value.  Mditionally, 
it  should  be  noted  that  the  initial  slope  of  the  response  is  not  zero  as  is 
true  for  a  second-order  system  without  a  zero  (13.2:91).  Ihe  block  diagram  of 
a  system  zero  is  shown  in  Figure  13.21. 


FIGURE  13.21.  BLDCK  DIfiGRAM  OF  A  SECOND-ORDER 
SYSTEM  WITH  A  REAL  ZERO 


Frcm  Figure  13.21  it  can  be  seen  that  the  zero  operates  on  the  input 
signal  to  produce  a  signal  proportional  to  both  the  magnitude  and  the 
derivative  (rate  of  change)  of  the  input  signal.  Therefore,  the  system  v/ill 
react  not  only  to  the  magnitude  of  input,  but  also  to  its  rate  of  change.  If 
R(s)  is  changing  r^idly,  then  E(s)  is  large  and  the  s^^^stem  responds  faster. 

(13.1:360) . 

13.7  STEADY-STATE  FREQUENCY  RESPONSE 

We  have  loo]:ed  at  the  time  domain  analysis  and  specifications  of  control 
systems.  In  the  time  domain  analysis,  the  ti'pical  test  inputs  were  the  step, 
ranp,  and  parabola.  The  frequency  response  technique,  introduced  in  this 
paragraph,  is  a  valuable  tool  to  the  control  systems  engineer  and  provides  a 
standardized  method  representing  the  total  performance  of  a  system.  The  input 
for  steac^-state  frequency  response  is  the  sinusoid 

r(t)  =  sin  u)t 

The  basis  for  the  frequency  response  method  is  that  a  system's  response 
to  a  sinusoid  will  be  a  sinusoid  at  the  same  frequency,  but  the  response  will 
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differ  in  magnitude  and  phase  angle.  All  that  is  needed  to  ccnpletely  specify 
the  steafy-state  frequency  response  is  to  be  able  to  find  the  magnitude  of  the 
output  and  the  phase  angle. 

Bie  feet  that  the  output  is  a  sinusoid  of  the  same  frequency  can  be  shown 
by  analyzing  a  sinusoid  iiput  to  a  first-order  system  described  by 

G(s)  =  (13.44) 

s  +  - 

T 


The  irput,  r(t)  =  sin  ojt  in  Laplace  transform  is 


R(s)  = 


(I) 


=,2.2 

S  +0) 


C(s)  =  G{s)  R(s) 


and 


C(s)  = 


B 


(1) 


12  2 
S  +  —  S  +  0) 

T 


Using  partial  feaction  expansion  yields 


C(s)  = 


S  +  -i 
T 


C2  S 

„2  .  2 
S  +  (1) 


=2.2 

S  +  0) 
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Finding  the  coefficients  C2,  and  can  be  a  tedious  process.  By 

inspection  can  write  the  form  of  the  solution  as 

-t/x  ^3 

c(t)  =  C,e  +  C-,  cos  tot  +  —  sin  ut 

1  Z  0) 

Another  form  of  this  equation  is 

c(t)  =  C^e  +  Aq  sin  (ujt  +  <|))  (13.45) 

Bie  stea^-state  response  can  be  written  as 


c(t)gg  =  Aq  sin  (ut  +  <|)) 

/  / 

^  Magnitude  * — 


(13.46) 


Phase  Angle 


viiich  tells  us  that  the  steady-state  response  will  always  have  the  same 
frequency  as  the  ii^jut  but  will  differ  in  phase  angle  and  magnitude.  Bie 

“t/  T 

transient  re^jonse  due  to  the  eaponential  term,  C^e  ,  decays  to  zero  as 
t-»«'. 

The  Laplacian  operator,  s,  contains  both  real  and  imaginary  conponents 
and  to  evaluate  coefficient  the  ccnplex  variable  "s"  would  be  selected  to 
be  (i.e.  purely  imaginary) .  Thus 

s  =  o  +  joj^ 

and  for  a  constant  aiiplitiade-irput  sinusoid,  a  is  zero,  (Figure  13.22) , 
therefore 


s  =  jto 


13.37 


Jw 


8-plane 


8  =  Ju) 

PURE  HARMONIC  MOTION 


(a  =  0) 


FIGURE  13.22.  S-PLANE  -  PURE  HARMONIC  MOTION 
The  frequency  response  function,  GCjco),  is  defined  by  replacing  s  with  jco 
in  the  system  transfer  function  (Equation  13.44). 

G(s)  «  — ^ 
s  +  i 


becomes 

G(jw)  -  — ^ 

j«  + 

It  is  inportant  to  remember  that  we  are  tallying  about  the  steady-state 
frequency  response  only  \i^en  we  replace  s  with  jw. 

13.7.1  Conplex  Numbers 

In  the  study  of  feedback  control  systems,  the  relative  magnitude  and  time 
relationship  between  such  quantities  as  position,  speed,  voltage,  current, 
force,  and  torque  are  the  items  of  interest.  These  are  all  real  physical 
quantities  v^ich  behave  according  to  the  laws  of  nature.  It  is  frequently 
convenient,  however,  to  represent  these  physical  quantities  by  ccmplex 
mathematical  symbols  that  indicate  more  than  the  information  describing  the 
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real  quantities  themselves.  The  use  of  conplex  verifies  to  represent  real 
fiiysical  quantities  has  the  advantage  of  sinplifying  the  mathematical  process 
necessary  to  solve  the  prvoblem.  On  the  othet  hand,  it  has  the  disadvantage 
of  obscuring  the  true  value  of  the  real  physical  quantities. 

It  is  the  purpose  of  this  chapter  to  introduce  the  conplex  variable 
notation  which  will  be  used  later.  Cooplex  quantities  are  usually  e3q)ressed 
in  one  of  four  forms: 

(a)  rectangular 

(b)  polar 

(c)  trigonometric 

(d)  exponential 

The  equivalence  of  these  four  forms  will  now  be  demonstrated. 

13.7.1.1  Rectangular  Form.  The  conplex  quantity  "z  is  drawn  on  the  conplex 
plane  in  Figure  13.23.  It  can  be  thought  of  as  a  position  vector  in  the 
conplex  plane. 


FIGURE  13.23.  THE  CCWPLEX  PLANE 

The  real  part  is  measured  along  the  horizontal  or  real  axis  and  the  imaginary 
part  is  measured  along  the  vertical  or  imaginary  axis. 


In  rectangular  form  the  complex  variable ,  z,  is 

7  =  X  +  jy 


where  j  is  the  imaginary  quantity  >/=r. 

13.7.1.2  Polar  Form.  Any  position  in  the  complex  plane  can  also  be  defined 
by  the  angle,  6,  of  the  position  vector  T,  and  its  magnitude,  |z’(( Figure 
13.23).  In  polar  form  the  complex  quantity,  'z,  is 


m 


In  terms  of  the  rectangular  form  parameters, 

iz|  .  v/7T7 

©  «=  tan”^  y/x 


13.7.1.3  Trigononaetric  Form.  The  trigonometric  form  of  the  position  vector 
in  the  complex  plane  can  be  written  again  using  Figure  13.23.  We  see  that 


cos  © 


and 


sin  © 


y 

\/?T7 


Therefore, 


cos  ©  +  j  sin  © 


Multiplying  both  sides  hy 


/ 


+  3 


X  +  y 


|2|  =  \/x^  +  y^ 
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#77  V?77j 

The  trigonanetric  form  is  then 

z  =  |z|  [cos  e  +  j  sin  6] 

13.7.1.4  Exponential  fOrm.  The  ejqxjnential  form  of  a  ccnplex  quantity  is 
most  convenient  for  mathematical  manipulation.  It  will  be  shown  equivalent  to 
the  trigoncrnetric  form. 

.  X  . 

.  The  MacLaurin  series  expansion  e  is 


|z|  [cos  6  +  j  sin  6] 


+  V 


=  1  +  X  + 


X 

21 


3^ 

3! 


+  . 


Letting  x  =  je 


02  03 

1  +  je  -  |r  -  j  It  + 


(13.47) 


Now,  sin  e  and  cos  6  can  be  defined  by  series  ejqsansions  as  follovre; 

3  5  7 

sine  =  (13.48) 

©2  6^  A® 

cose  =  (13.49) 


Recalling  that  j  =  -1,  Equation  13.48  and  13.49  be  written 


oos  e  =  1  +  iMi  +  JmI  +  iiltl  + 

cot.  o  i-r  91^  41 


(13.50) 
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j  sin  e 


(13.51) 


Mding  Equation  13.50  to  13.51  yields 


cos  6  +  j  sin  e  =  1  +  4t-  +  . +  +  .  . 


li 


21 


3! 


(13.52) 


Bie  right  side  of  Equation  13.52  is  equal  to  Equation  13.47,  therefore 

e^®  =  cos  e  +  j  sin  6 


and  finally 


2  =  |z|  e 

We  have  proven  the  four  forms  of  the  coiplex  variable  z 
Biey  are  sunmarized  below 


(13.53) 
to  be  consistent. 


Rectangular 

Polar 

Trigononetric 

E35»nential 


z 


z 


z 


z 


X  +  jy 


[cos  6  +  j  sin  6] 
e^® 


13.7.2  Bode  Plotting  Itechnique 

Willi  this  bacl^ground  in  ccnplex  notation  we  will  develop  the  Bode 
technique  of  frequency  response.  Beginning  with  a  generalized  transfer 
function  (Equation  13.54),  we  will  manipulate  it  into  the  frequency  response 
standard  form,  sometimes  called  the  Bode  form  (Equation  13.55) . 
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G(s) 


K  s*  (s  +  aj(s  +  a  )(s^  +  2C(o 


)„  S  +  (0  ) 

n ,  n  ' 

1  1 


s  (s  +  aj)(s  +  a.  )(s^  +  IZfn  s  +  o)  * ) 


(13.54) 


13.54  must  be  normalized  as  follows: 


G(s) 


~  ^  1  \ 

_  ^  j 


^3  ^4  «  "  s" 


(•^  s  +  1)  (t.  S  +  1)/^  -+ 


<*L 


(13.55) 


idiere 


Let 


I  »  etc, 
“2 


®3  a, 


and  K  -  static  loop 
sensitivity 


SubsUtute  jM  for  s  in  Bjuation  13.55  and  rewriting  in  the  Freijiency  Besponse 
Standard  Form 


K  (j«)"  (1 


G(ja)) 


+  jTjW)  (1  +  jTjW)  I  1  -  -id  +  j22;  ^ 


( j«)"  (1  +  jTjW)  (1  +  jT  «)  /  1  - 


~  +  j2C 
0) 


0) 

«n 
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Equation  13.56  can  be  written  in  polar  form 

G(jco)  =  |G{j»)| 

v^ere  jG(ja>)  j  is  of  the  form  +  Im^ 

and  ♦(«)  of  the  form  tan  ^ 

Furthermore,  G(  jw)  can  be  written  in  eaqponentied  form  where 

G(jw)  «  |G(jw)|  (13.57) 

TO  egress  G(j»)  in  either  of  these  formats  will  entail  finding  the  magnitude 
G(j«)  and  the  phase  angle  ♦((•)).  The  Bode  technique  requires  taking  the  log  of 
G(j<i))  to  take  advantage  of  addition  and  subtraction  in  lieu  of  multiplication, 
and  division. 

Taking  the  log  of  Equation  13.57  yields 

log  |G(ja))|e^*<“>  =  log  |G(ja))|  +  log 

-  log  |G(j«)|  +  j  .4343  ♦(«) 

The  qucuitity,  j  .434  ♦(&)),  is  the  imaginary  part  and  in  future  discussion 
only  the  angle  ^((o)  will  be  used. 

In  feedback  system  work  the  unit  conDnonly  used  for  the  logarithm  of  the 
magni tilde  is  the  decibel.  The  logarithm  of  the  magnitude  of  a  tremsfer 
function  G(ja))  expressed  in  decibels  is 

20  log  |G(jfi))|  =  (  )  db  (13.58) 

This  quanti^  is  often  referred  to  as  the  log  magnitixie  and  is 
abdsreviated  Lm  v^ere 


Lm  G{  j(o)  =  20  log  |G(j<o)|  db  (13.59) 
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Now,  how  does  multiplication  and  division  become  addition  and  subtraction 
for  Bode  development?  We  will  take  the  log  of  Equation  13.56  and  multiply 
this  by  20  \^ich  will  give  the  anplitude  ratio  in  decibels.  Bie  use  of 
logarithms  will  edlow  us  to  add  for  multiplication  and  subtract  for  division. 

The  Lm  of  Equation  13.56  becomes 


20  log  |G(  j»)  I  -  20  log  +  20  m  log  |  j«|  +  20  log  |l  +  j  Tjwj 
+  20  log  |1  +  jTjWl  +  20  log  |1 +  j  | 

“a. 


-  20  n  log  |ja)|  -  20  log  jl  +  jT-wj 


;  2Cft> 


-  20  log  |1  +  jT^«|  -20  log  l--~  +^1^ 


(13.60) 


The  associated  phase  angle  of  Equation  13.56  becomes. 


+  m  90®  +  tah"^  oxt^  +  tan'^wTj 


2Cco/a> 


+  tan“ 


1  - 


a> 


-  n  90®  -  tan”^ 


A) 

n 


2CA>/b) 


-  tan“^  (OT^-  tan~^ 


1  - 


CO 


CO 
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It  is  iirmediately  cibvious  that  a  point'  by  point  solution  of  these 
equations  for  varying  ii^ut  frequency,  w,  would  be  very  tedious.  A  technique 
known  as  the  Bode  plot  simplifies  this  process.  Notice  in  Equation  13.60  that 
there  are  four  types  of  factors  in  the  open— loop. transfer  function  G(ja)) . 


1.  Constant  term,  K 

n 

2.  Pole  or  zero  at  origin,  (ju) 
(+n  =  zero,  -n  =  pole) 


±n 


+n 


3.  Siitple  pole  or  zero,  (1  +  jux)  — 

/  2  (D 

4.  Quadratic  pole  or  zero,  j  1  -  +  j2c  ~  I 

V  0,^  n/ 


Ihe  Bode  plot  uses  semilog  paper.  Magnitude  and  piiase  angle  are 
represented  on  the  ordinate  (linear  scale)  and  frequency  along  the  logarithmic 
scale  as  in  figures  13.24  and  13.25.  Bode  plot  technique  uses  asyirptotes  and 
corrections  to  the  a^ptotes.  for  each  of  the  four  types  of  factors  listed 
above.  All  of  the  factors  are  individually  plotted  on  the  Bode  diagram,  and 
then  are  and  subtracted  (taking  advantage  of  logarithms)  to  achieve  the 

conposite  curve.  We  will  develop  the  technique  for  each  of  the  four  types  of 
factors. 

Constant  term,  Ihe  magnitude  of  in  db  is 

20  leg  j  j  =  constant 
and  the  associated  phase  angle  is 

Arg  (K^)  =  0°  or 

Arg  (-  K^)  =  +180° 

as  shown  in  Figure  13.26. 

Bie  magnitude  and  jiiase  angle  are  depicted  respectively  on  Figures  13.24 
and  13.25. 
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PHASE  ANGLE  -  DEGREES  MAGNITUDE  -  DECIBELS  (db) 


FIGURE  13.26.  PLOT  OF  A  CONSTANT  ON 
<niE  COMPLEX  PLANE 

Pole  or  zero  at  origin,  { jw)  aagnitude  of  ( jw)  *  is 

20  log  |(ja))  |  -=  ±20  n  log 

is  the  equation  of  a  straight  line  with  sl(^  of  ±20  n  db/decade. 

A  decade  is  a  frequency  band  from  fj  to  fj  idiere  fj/fj  “  10*  If  fj  ■  4 

rad/s,  then  one  decade  above  f^  is  :  f^  -  40  rad/s.  The  octave  is  also  used 
as  a  frequency  ratio  and  is  a  frequency  band  from  f^  to  fj  %diere  fj/f^  «  2. 
If  ■  4  ta6/s,  then  one  octave  above  f^  is  :  fj  =  8  rad/s.  The  slope  of 
+20  n  db/decade  -  ±6  n  db/octave. 

When  (0  »  1»  the  term  20  log  j  (jw)  *"|  *  20  log(l)  »  0.  Therefore,  the 
magnitude  plot  of  a  pole  or  zero  at  the  origin  is  a  straight  line  with  a 
slope  ■  ±  20n  dB/decade,  passing  through  the  0  db  point  at  frequency 
ci>  =  1  (Figure  13.24).  The  ptese  angle  of  (jw)  *"  is 

Arg|(jft))  «  ±n  90® 


as  shown  in  Figure  13.25. 


Simple  Pole  or  Zero,  (1  +  jcor)*”  the  magnitude  of  this  term  in  db  is  expressed 
as 


20  log  j(l  +  jwT)  1  «  ±20  n  log  (1  +  j«t) 
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LOG  MAGNITUDE  -  DECIBELS  (db) 


at  very  lav  fraguency  (i.e.,  wt  «  1)  the  magnitude  of  this  curve  is  0  db.  At 
fe-eguencies  vAiere  wt  »  1  the  magnitude  asynptote  has  a  slc^  of  +20n 
db/decade.  Ihe  0  db'  asymptote  and  the  +20n  db/dec  asymptote  intersect  at  the 
corner  frequency,  =  l/x,  Figure  13.27. 

Ihe  phase  angle  is  sqrressed  as 

Arg  (1  +  jwT)  ^  =  ±  n  tan”^  ut 
At  0)  =  0,  <|)  =  0°  and  at 

=  ±n  t/2  radians 


FIGUPE  13.27.  BODE  !©GNITUDE  PLOT  OF  TEEM  (1  +  jtox)  - 


Table  13.2  shows  the  variaUon  of  the  phase  angle  with  nonnalized  frequency 
w/Wp,  for  n  =  -1. 

TKBI£  13.2 

PHASE  ANGLE  VAKEATION  WITH  NOIWALIZED  FREQUENCE  Hi  +  jwT)"^! 


(0 

"c 

tan”^  uT 

0 

0 

.1 

-5.7° 

.5 

-26.5° 

1.0 

-45.0° 

2.0 

-63.4° 

10.0 

-84.3° 

GO 

-90.0° 

in  » 

The  following  techniques  are  used  to  plot  the  (1  +  jwx)  factor. 


1. 

2. 


Locate  comer  freguency,  o)^  -  I/t. 

Draw  +20n  db/decade  asynptotes  through  the  comer  frequency  (+20n 
db/dec  for  zero  temis  and  -20n  db/dec  for  pole  terms) . 

A  straight  line  can  be  used  to  apprcodmate  the  phase  shift.  The 
line  is  drawn  fron  0°  at  one  decade  below  the  com^  f requai^  to  n 
(+90°)  (+  for  zero  term,  -  for  pole  term)  at  one  decade  ah^e  th 
comer  frequency.  The  maximum  deviation  using  this  approso^tron  is 
about  6°.  The  specific  phase  angle  v^ues  in 

and  the  ^propriate  corrections  can  be  ^>plied  if  desired.  Tne  e 
cx)rrections  are  shewn  in  Figure  13*28. 

The  error  to  the  magnitude  curve  (created  by  using  the  asymptote 
technique)  can  be  determined  analytically.  First  determine  the 
error  at  the  comer  frequency  —  l/i. 


+20n  log 


1  +  0)  ^  becones  +20n  log 
c 
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and 


+20n  log  vr  =  +10n  log  2  =  +3.01n  db 

Biis  shows  that  the  asymptote  can  be  corrected  adding 

aFproKiinately  +3n  db  at  the  comer  frequencjr.  Likewise  for  a 
frequency  one  decade  above  the  comer  frequency. 


(I)  =  10  w  ~ 

C  T 


and 


+20n  log 


^f 


1  +  100 


2/  2 
T  /t 


+20n  log  VloT 


and 


-10  log  101  =  -20.043  db  (actual,  for  n  -  -1) 

Our  straight  line  asynptote  used  -20  db  so  the  total  error  at  one  decade  is 
-.043  db.  Similarly  the  error  at  w  /lO  can  be  found.  At  w  =  2a)  (one 
octave)  =  2/t,  the  actual  Im  for  n  =  -1  is 

-20  log  -^1  +  4  =  -10  log  5  =  -6.9897  db 

The  asyiiptote  method  produced  a  value  of  -6  db  at  this  point,  thus  an  error 
-.9897,  or  appraxiraately  -1.0  db. 

Therefore,  the  strai^t  line  asymptotes  can  be  made  closer  to  the  actual 

Im  curve  by  applying  a  +3n  db  correction  at  o)^  and  a  +ln  db  correction  1 

octave  above  and  1  octave  below  o)  . 

c 
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FIGURE  13.28.  BC»E  PHASE  ANGLE  PLOT  OF  (1  +  joit)  ^ 
Qaadratic  Tterniy  (1  +^c/a)^  ju  + 

Oonsider  the  quadratic  term 


G(s) 


1 


+ 


(13.61) 


Hie  log  magnitude  of  G(ja)),  Equation  13.61,  in  db  is 


20  log  G(ju)) 


(13.62) 


and  tile  phase  angle  is 
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FIGURE  13.29.  BODE  DIAGRAM  FOR  G(j(i)) 


1  +  jo)  +  (ju/Wj^) 
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At  very  low  frequency  ^  «  1 


20  log  G(j(i>)  s  -  20  n  log  1  «  o  db 

Therefore,  at  low  frequencies  the  asymptote  for  the  quadratic  term  is  a 
straight  line  with  a  slcpe  of  zero.  At  high  frequency  »  1,  Equation 

13.62  can  be  approximated  as 


20  log  G(j(o)  «  -20  n  log 


as  -40  n  log  ^  ^  ^  db 


The  last  equation  represents  the  equation  of  a  straight  line  with  slope 
of  -40n  db/dec.  If  the  quadratic  is  in  the  numerator  (i.e.,  +n)  the  slqpe  is 
positive.  The  two  asymptotes  intersect  at  (o  ,  hence,  o)  is  considered  the 
comer  frequency  of  the  quadratic  factor.  The  actucd.  magnitude  plot  for  the 
quadratic  factor  differs  strikingly  from  its  asymptotic  plot  in  that  the 
amplitude  curve  depends  not  only  on  but  also  on  C.  Prom  Figure  13.29, 
several  values  of  im  around  can  be  plotted  for  a  specific  C  to  obtain,  an 
accurate  magnitude  plot. 

The  phase  angle  plot  for  a  guadratic  factor  can  be  obtained  locating 
the  +90 ®n  point  at  and  obtaining  a  few  points  either  side  of  for  a 
specific  value  of  C  frcmn  Figure  13.29. 
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To  suntnarize,  the  procedures  for  plotting  the  queidratic  term  are: 


a.  Determine  the  value  of  c  and 

b.  Plot  the  zero  db  asymptote  fron  lc3w  frequencies  to  and  a  +40n 
db/dec  asynptote  beginning  at 

c.  Use  the  curves  presented  in  Figure  13.29  to  correct  the  eisymptotes 
in  the  vicinity  of  the  comer  frequency. 

d.  At  the  comer  frequency,  locate  the  ±90*^n  phase  point.  Using 

the  cujnres  of  figure  13.29  for  the  specific  c,  plot  enough  data 
points  to  permit  sketching  the  phase  angle  curve. 

When  each  of  the  four  types  of  factors  are  plotted  on  tiie  Bode  plot,  all 
the  magnitude  curves  and  phase  angle  curves  are  summed  at  different 
■Fraty^^'nci  to  Complete  the  composite  curves.  Ihe  follcMing  problem  vdJ.1 
illustrate  the  simplicity  of  this  technique. 

640s  (s  +  1000) 

(s  +  10) (s^  +  80s  +  6400) 

2cu^  =  80,  03^  =  80,  C  =  0.5 
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2.  Find  the  cx>mer  frequencies  where  a)_  ^  l/x.  For  a  quadratic  term 
the  natural  frequency  is  the  comer  frequency 

Zeros:  u  =  1000,  +20  db/dec 

C 

Poles:  (0  =  10,  -20  db/dec;  w  =  80,  -40  db/dec 

c  c 

3.  Plot  the  individual  magnitude  and  phase  angle  terms  on  the  Bode. 
Also  20  log  =  20  log  10  =  20  db. 

4.  i^ly  the  appropriate  corrections  at  the  comer  frequencies. 

5.  Add  the  curves. 

6.  Figures  13.30  and  13.31  sha*;  the  contribution  of  the  separate 
&cbors  and  the  oonposite  curve  for  the  example. 

We  must  eniphasize  that  the  development  on  the  Bode  plot  presented  here  is 
based  on  the  stea^-state  frequency  response  of  an  open-loop  system  to  a 
sinusoidal  input.  Techniques  exist  to  arrive  at  the  closed-loop  feequency 
re^x>nse,  but  these  are  beyond  the  scope  of  our  stu<^.  Ihe  closed-loop 
frequency  response  graph  is  a  plot  of  magnitude  ratio,  M(ja))  =  C(ja))/R(jo)) , 
and  phase  angle,  <|»,  versus  frequency.  One  method  of  determining  the 
closed-loop  frequency  response  is  by  tising  the  Nichol's  chart  (Reference 
13.3). 
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FIGURE  13.30.  BODE  LOG  MAGNITUDE  PLOT 


FIGURE  13.31.  BODE  PHASE  ANGLE  PLOT 


13.7.3  Relative  Stability 

file  relative  stability  of  a  closed-locp  system  can  be  determined  by 
looking  at  the  Bode  plot  of  the  open-loop  transfer  function,  KGBKjw) .  Several 
terms  are  vised  to  relate  stability  by  the  Bode  Plot.  The  mathematical  basis 
for  these  relationships  comes  from  the  Nyguist  Stahili-ty  Criteria.  The  terms 
are: 

13.7.3.1  Gain  Margin.  Gain  margin  is  the  additional  amount  of  gain,  measured 
in  decibels,  that  the  magnitude  ratio  can  be  increased  before  the  system  goes 
unstable.  The  gain  margin  is  defined  as  the  reciprocal  of  the  open-loop 
transfer  ftmction,  GHCjco),  evaluated  at  the  firequenoy  vbere  the  phase  angle  is 
-180°. 


Gain  Maigln  =  20 


This  quantity  is  illustrated  in  Pigure  13.32. 

13.7.3.2  Phase  Margin.  Phase  margin  is  the  amount  of  phase  shift,  measured 
in  degr^s,  that  the  phase  angle  curve  can  be  displaced  to  produce  instability 
in  the  system.  Phase  margin  is  measured  at  the  frequency  vbere  the  Im  plot 
crosses  the  0  db  line. 

Phase  Margin  =  +180°  +  «|) 

(4i  =  phase  angle  measured  at  0  db) 

This  quantity  is  illustrated  in  Pigure  13.32. 
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stability  requires  that  the  ^Aiase  margin  be  positive,  i.e. ,  the  phase  angle  at 
the  0  db  crossover  point  must  be  greater  than  -180^. 

13.7.3.3  Bie  Gain  Crossover  Point  is  defined  as  the  point  or  points  vdiere  the 
magnitude  curve  crosses  the  zero  db  axis. 

13.7.3.4  The  Phase  Crossover  Point  is  the  point  on  the  Bode  phase  angle  plot 
at  vdiich  the  phase  angle  is  -180°.  Ihe  frequency  at  vMch  the  phase  crossover 
occurs  is  called  the  gain  nazgin  frequency. 

13.7.4  Frequency  Domain  Specifications 

%ere  are  several  tezms  to  express  the  specifications  of  systems  in 
the  frequency  domain.  Althou^  these  tezms  are  usually  used  to  define  the 
closed-loc^  response,  they  can  also  be  used  to  express  characteristics  of  the 
open-loop  Bode  plot. 

13.7.4.1  Bandwidth  (BW) .  The  definition  of  bandwidth  of  a  ^tem  depends  on 
an  accurate  description  of  the  problem.  Nozmally  the  bandwidth  is  defined  as 
the  frequency  at  \diich  the  magnitude  ratio  M(ja))  =  C(jo))/R(jto)  has  dropped 
to  70.7%  of  the  zero  frequency  level  or  3  db  down  from  the  zero  frequency 
level  as  shown  in  Figure  13.33.  This  does  not  cover  all  cases  in  that  the 
magnitxide  ratio  at  zero  frequency  may  be  low  as  in  Figure  13.34.  In  this  case 
the  bandwidth  is  defined  as  the  frequency  range  over  vdiich  the  magnitude  ratio 
does  not  vary  more  than  -3  db  frcm  its  value  at  a  specified  frequency.  For  the 
purposes  of  this  ccxirse  the  bandwidth  can  also  be  determined  as  the  -3  db 
point  on  the  Im  plot  of  the  open-loop  systan.  This  value  should  correspond 
closely  to  the  bandwidth  of  the  closed-loop  system. 
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|M|  -  dbs 


BANDWIDTH 


Odb 

-3db 


^  I 

_ _ I  I 

«r  "co 

FREQUENCY 


FIGURE  13.33.  FREQUENCY  DOMAIN  CHARACTERISTICS 


FIGURE  13.34.  FREQUENCY  DCmiN  CHARACTERISTICS 
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The  frequency  at  which  the  -3  db  point  is  reached  is  called  the  cutoff 
frequency/  ^co*  Bandwidth  is  uiportant  for  two '  reasons.  First/  it  is 
indicative  of  the  noise  filtering  characteristics  of  the  ^tem.  Syston  noise 
is  always  present  and  the  bandwidth  and  the  corresponding  cutoff  frequenq^ 
dictate  at  what  frequency  the  response  and  thus  the  noise  will  be  filtered. 
Secondly/  the  bandwidth  is  a  raeasure  of  the  transient  response  properties  of  a 
Qfstan.  A  large  bandwidth  will  allow  hi^er  frequencies  to  pass  to  the  output 
and  the  system  may  be  characterized  by  faist  rise  time  and  large  overshoots. 
However/  if  the  bancb^idth  is  narrow,  only  low  frequency  signals  are  passed  and 
the  response  will  generally  be  slew  and  sluggish. 

13.7.4.2  Itesonant  Ffeak,  M^.  If  the  ^stem  is  of  seoond-order  or  hi^er,  it 
may  have  a  resonant  peak,  M^.  For  a  secjond-order  system  there  exists  exact 
mathematical  relationships  betveen  5,  the  daiiping  ratio,  and  the  frequent 
at  v^iich  M  occurs.  A  hic^er-order  systen  can  often  be  aH>rcximated  by  a 
second-order  ^stem  to  simplify  the  solution.  Ihe  resonant  peak,  Figure 
13.33/  is  an  indication  of  the  relative  stability  of  the  system  as  a  high 
value  of  <r>rr«=>sp»iTdfi  to  a  large  overshoot  in  the  time  demain.  lypical 
values  of  for  usable  stable  ^sterns  may  ■vary  from  1.1  to  1.5. 

13.7.5  Bqjerimental  Method  of  Frequency  Response 

A  Bode  plot  m^  be  determined  experimentally  and  \iLtiroately  will  provide 
"the  system  transfer  function.  Ihe  method  depicted  in  Figure  13.35  will  allow 
for  measurement  of  the  magnitude  ratio  and  phase  angle  versus  frequency. 


SINE  WAVE 


BRUSH  RECORDER 


FIGURE  13.35.  EXPERIMENTAL  BCDE  TECHNIQUE 
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Once  the  Lm  and  phase  angle  versus  frequency  curves  are  plotted,  asymptotes  can 
be  fitted  on  the  curve  and  the  corner  frequencies  determined.  If  a  resonant 
peak  occurs,  use  the  techniques  discussed  previously  to  determine  2;  and  («>  . 

n 

In  our  discussion,  we  have  been  talking  about  open-loop  systems  that  do 
not  have  poles  and/or  zeros  in  the  right-half  of  the  s-plane  (RHP).  These 
systems  are  known  as  minimum  phase  systems.  A  nonminimum  phase  system  is  one 

which  has  an  open-loop  pole  and/or  zero  in  the  RHP.  A  nonminimum  phase  factor 
is  of  the  form 


(l-jC»)*"  or  1+^2^ 

“n  ■ 

If  a  system  is  known  to  be  minimum  phase,  only  the  Lm  plot  is  required  to 

fully  determine  the  i^stem  transfer  function;  while  both  the  Lm  and  pdiase 

angle  plots  are  required  to  determine  the  system  transfer  function  if  the 

^tem  is  nonminimum  phase.  For  exanple,  consider  the  following  transfer 
functions: 


Gj  (s) 


(s+1)  , 

TiOTT 


G2(s) 


(s+1) 


'  G3(s) 


(s-1) 

(s+lO) 


'  and  (s) 


(s-1) 

Ti^ 


All  four  transfer  functions  have  similar  Lm  plots,  but  their  phase  angle  plots 

are  all  different  (Gj(s)  is  minimum  phase;  G2(s),  G3(s),  and  G4(s)  are  all 
nonminimum  phase). 


13.8  CLO5EE>-L0OP  TRANSFER  FUNCTICN 


For  reasons  to  be  seen  shortly,  conplex  control  systems  are  most  often 
represented  by  a  block  diagram  in  the  form  of  Figure  13.36  \diere  the  forward 
transfer  function,  G,  and  the  feedback  transfer  function,  H,  are  expressed  as 
functions  of  s,  the  Laplace  transform  variable.  The  closed-loop  transfer 
function  of  Figure  13.36  will  now  be  developed  in  terms  of  the  forward  and 
feedback  transfer  functions  for  our  first  and  second-order  systems.  No 
transfer  function  present  in  the  feedback  loop  is  called  unity  feedback  (in 
Figure  13.5,  H(s)  =  1). 
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FORWARD  TRANSFER  FUNCTION  —  G(s) 

FEEDBACK  TRANSFER  FUNCTION  -  H(s) 

OPEN-LOOP  TRANSFER  FUNCTION  (OLTF)  —  G(s)  H(s) 
CLOSED-LOOP  TRANSFER  FUNCTION  (CLTF)  -  C(s)/R(s) 
ACTUATING  SIGNAL  —  E(s) 


FIGURE  13.36.  STANDARD  FORM  CF  FEEDBACK  OCXJTROL  SYSTEM 

In  order  to  find  G(s)  of  the  first-order  ^tem  (Figure  13.5) ,  equate 
torques  and  assume  only  danping  present 

ve  =  b  6 

taking  the  Laplaoe  transform  and  noting  that  G(s)  =  e(s)/E(s)  we  find 

yE(s)  =  bse(s) 

G(s)  =  I  (s)  =  ^ 
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again  letting 


*  4  (13.66) 

Uie  sane  procedure  for  the  case  including  inertia  yields  the  follovang 
torque  summation: 


fjz  ■=  10  +  b  6 


E(s)  »  ^  e(s)  +  ^  e(s) 


G(s)  -  I  (s)  -  —IL — 

is^  +  bs 


(13.67) 


Thus^  vre  have  the  forward  transfer  function  for  our  two  systems.  Referring  to 
Figure  13.36,  we  will  now  derive  an  e:q)ression  for  the  closed-loop  transfer 
function  in  terms  of  G  and  H. 


G(s) 


C(s) 

E(¥T 


and  also 


E(s)  =  R(s)  -  H(S)  C(S) 
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substituting 


G(s)  [R  -  H(s)  C(s)]  =  C(s) 


R(s)  G(s)  =  C(s)  U+  GH(s)]  (13.68) 


the  system  closed-loop  transfer  function  beccmes 


G(s) 

1  +  GH(S) 


(13.69) 


This  is  a  very  iirportant  relationship  which  should  intnediately  be  ccnmitted  to 
memory. 

As  noted  previously,  the  block  diagram  of  Figure  13.36  is  the  standard 
form  of  the  feedback  control  ^stem.  When  in  this  form  the  closed— loop 
transfer  function  can  be  quickly  found  by  Equation  13.69.  But  most  irtportant, 
the  characteristic  equation  of  the  system  from  v^ch  the  transi^t  response  is 
determined  is  inmediately  evident.  Referring  to  Equations  13.68  and  13.69  we 
will  show  that  the  characteristic  equation  is  found  from  the  denominator  of 
the  ri^t  hand  term  in  Equation  13.69. 

1  +  GH(s)  =  0 
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The  charac±eristLc  eguald.on  is  merely  1  plxis  tiie  system  open-loop  transfer 
function,  GH(s) ,  vjhich  is  directly  available. 

But  first,  applying  our  closed-loop  transfer  function  expression  to  the 
first-order  i^stesm  we  rase  the  forward  transfer  function  of  Equation  13.66. 
Since  the  system  has  amit  feedback,  H(s)  =1  and 


GH(s) 


TS 


Therefore,  using  Equation  13.69  yields 


2^ 

£  =  G(s)  _  TS 

R  1  +  GH(s)  “  1  +  1_ 

TS 


1 

TS  +  1 


vdiich  is  consistent  with  Equation  13.16,  the  transfer  function  derived  fmm  the 
equation  of  motion  of  the  entire  system. 

In  the  case  of  the  second-order  ^stero,  since  H{s)  is  uni1^ 

GH(s)  =  — -H - 

Is^  +  bs 


For  this  ^^stem,  axsing  Equation  13.69 


_ y 

£  _  G(s)  _  Is^  +  bs 

R  1  +  GH(s)  "TT  u 

2 

Is^  +  bs 


C 

R 


(s) 


1 


1 


(13.70) 


is  3.93.in  consistent  with  the  more  direct  method  leeding  to  Equation  13 .IS. 
Ihe  dencminator  of  Equation  13.70  is  the  characteristic  equation  introduced  by 
Equation  13.12. 

Transfer  functions  are  written  to  describe  either  vAiole  systems  or  parts 
of  systems  usinq  the  appropriate  differential  equation.  When  control  systems 
described  by  block  diagrams,  are  reduced  to  some  standard  form,  they  quickly 
yield  both  the  transfer  function  of  the  entire  closed-loop  system  and  its 
characteristic  equation. 

Wte  will  now  discuss  the  technique  of  manipulating  control  systems  in 
block  diagram  notation  to  obtain  the  desired  form. 

13.9  BLOCK  DIftSRAM  ALGEBRA 

It  was  seen  that  the  sirtplificaticai  resulting  fron  the  use  of  operational 
cal  ml  ns  is  further  increased  vdien  transfer  functions  and  block  diagrams  are 
introduced.  Ihe  special  methods  of  predicting  the  transient  response  of  a 
system  withcut  solving  its  equation  of  motion  are  most  conveniently  employed 
vAien  the  block  diagram  is  of  the  form  of  Figure  13.36. 

In  practice,  individual  transfer  functions  are  written  for  each  integral 
unit  of  a  more  ocirplex  ^tem.  For  example,  the  system  of  Figure  13.37 
represents  the  pitch  axis  of  an  aircraft  autopilot  vhere  the  input  is  the 
ccromanded  pitch  attitude  and  the  outqait  the  actual  aircraft  attitude.  The 
autopilot,  the  elevator  servo,  and  the  aircraft  itself  are  described  separately 
in  G2,  and  G^  respectively.  As  long  as  it  is  realized  that  transformed 
quantities  are  used  the  G(s)  can  be  discarded  and  only  G  xjsed. 

The  system  of  Figure  13.37  can  be  simplified  by  coiibining  the  inner  loop 
into  a  single  transfer  function.  If  we  let  G^  be  the  closed-loop  transfer 
function  of  the  inner  loop  we  have 

S  “  1  +  G2 

Figure  13.37  can  then  be  redrawn  as  shown  in  Figure  13.38.  This  diagram  is 
then  further  reduced  by  noting  that 
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Dsnoting 


6.6  P 

«i  =  r-  =4  =  67  '  =3  =  C 

X  e 


6.  6  - 

rrr  - 

E  6.6  E 
i  e 


s  =  Ws  =  i¥4 


AIRCRAFT 

AUTOPILOT  DYNAMICS 

V  ELEVATOR  SERVO  j 


7 


R  +/OvE  I  G,  I  ^ 


•  Gj  - 63 


ACTUATING  /  \ 

\ 

SIGNAL  /  \ 

>- ELEVATOR 

I  L  N 

ANALF 

^COMMANDED 

^  ELEVATOR 

ELEVATOR 

ACTUATING 

ANGLE 

SIGNAL 

FIGURE  13.37.  AIRCRAFT  PITCH  AXIS  OCNTROL  SYSTEM 
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FIGDEE  13.39.  (FIGURE  13.37  FURTOER  REDUCED) 

We  have,  finally,  the  cxintrol  system  described  in  the  proper  form  in  Figure 
13.39. 

Hie  following  block  diagram  identities.  Figure  13.40,  will  assist  in 
manipulating  cxirplex  control  systems  into  the  standard  form  for  analysis. 
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4 


FIGURE  13.40.  CCRT.  BIDC3(  DIA3RAM  IDEKTITIES 


13.10  STEADY-STATE  PERFCMMaNCE 

The  stea^-state  accuracy  of  a  system  is  of  considerable  irtportance  and 
is  often  related  in  terms  of  the  stea%-state  error,  figures  of  merit  for 
stea<fy-state  performance  are  the  error  constants,  K_,  K  ,  and  K  often 
referred  to  as  the  position,  velocity  and  acceleration  error  constants. 

A  technique  used  to  indicate  the  steac^'-state  performance  of  a  i^stem  is 

to  classify  the  system  by  "Type".  The  number  of  free  or  pure  integrators  in 

the  forvrard  loop  is  the  Type  ^stem  (systan  must  be  stable  and  represented  by 

imity  feedback) .  For  a  specified  iiput  function,  an  n-type  system  will 

produce  a  mathematically  predictable  steacty-state  error.  Consider  the  unity 

feedback  system  in  Figure  13.41.  A  unity  feedback  system  is  lased  in  this 

development  since  we  will  be  relating  the  performance  as  a  function  of  the 

steac^-state  error,  e(t)  ,  v^ere 

ss 


e(t)^  - 


For  this  relationship  to  be  valid  in  this  development,  the  reference  iiput 
r(t)  and  the  control  variable  c(t)  must  be  dimensionally  the  same  and  must  be 
to  the  same  scale. 
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FIGURE  13.41.  UNITY  FEEDBACK  SYSTEM 


E(s)  =  R(s)  -  C(s) 


and 


C(s)  =  E(s)  G(s) 


therefore 

E(s)  =  R(s)  -  E(s)  G(s) 

E(s)  [1  +  G(s)]  =  R(s) 

E{s)  _  1 

R(s)  1  +  G(s) 

The  error  signal  E(s)  is  a  function  of  the  plant,  G(s) ,  and  the  reference 
input  R(s) . 

G(s)  can  be  represented  by 


G(s)  = 


K„(T  s  +  1)  (t  S  +  1)  -- 
n  1  1 


n  / 


(T  S+  1)(T  =+  1 

I  u  2 


n  cm 


■) 


Defines  "Type" 
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viiere  n  =  0,  1,  2,  ....  Ftor  a  0  system  n  =  0;  i.e. ,  no  free 

integrators  in  the  forward  loop.  G(s)  must  be  e:^ressed  in  the  above  form  to 
properly  evaluate  the  overall  gain,  K^,  of  the  transfer  function.  Uiis  gain 
is  often  referred  to  as  the  "DC"  gain  or  "l^^pe"  gain. 

We  will  se^  to  show  the  relationships  between  the  inputs,  the  n-type 
system  and  the  stea<fy-state  error,  first  cxDnsider  the  general  error  and  ^^ly 
a  step,  ranp,  and  a  parabolic  ii^t. 


E(s) 


R(s) 

1  +  G(s) 


13.10.1  Step  Input 

Let  r(t)  =  Ru_j^(t)  and  R(s)  =  R/s  and  apply  the  final  value  theoran. 
Recall  for  any  ftmction  F(s) . 


lim  s  F(s)  =  lim  f(t)  E(s) 

s  0  t  » 


lim 

s-»-0 


s  R/s 
1  +  G(s) 


R 

1  +  lim  G(s) 
s^O 


1 

1  +  G(s) 


R(s) 


K  Position  Error  Oonstant.  Ihe  position  error  constant,  I^,  is  defined  as 

'  '  ' 


K  =  lim  G(s) 
^  s-»-0 


Therefore, 


e(t) 


ss 


R 
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Far  a  Type  0  system 


G(s)  =  K, 


( - )  ( - ) 

0  ( - )( — ) 


K_ 


lim  G(s) 
s-»-0 


=  lim  K, 
s-*-0 


( - )  ( ) 

0  ( - )( ) 


K  ~  overall  gain  of  the  transfer  function. 


Ftor  a  Type  0  system  a  step  input  yields 


e(t) 


ss  1  +  K, 


and  the  steacty-state  error  is  represented  grajMcally  in  Figure  13.42 


(STEP  INPUT) 


- /- 


\y 


T 


c(t) 


;  (STEADY 
/  STATE  RESPONSE). 


1  +Ko 


TRANSIENT 


FIGURE  13.42.  STEADY-STATE  ERROR, 

TYPE  0  SYSTEM  -  STEP  INPUT 


For  a  Type  1  system 


G(s)  — 


K^( - )  ( - ) 

S  ( - )( - ) 
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K  =  lim  G(s) 
^  s-»-0 


=  lim 
s^O 


Ti=n=r 


K 


Therefore, 


e(t) 


ss  1  + 


=  0 


The  resulting  error  from  a  step  ii^jut  into  a  Type  1  system  is  zero.  Similarly 
for  a  Type  2  and  3  system  with  a  step  input,  the  resiilting  error  is  zero. 

13.10.2  Ranp  Input 

2 

Consider  a  ranp  iiput  r(t)  ==  Ru_j^(t) ,  R(s)  =  R/s 
Therefore, 


E(s)  =  ^ 


1  +  G(s) 


s-»-0 


s  R/s^ 


1  +  G(S) 


=  lim 
s-»-0 


R 


s  +  s  G(s) 


lim  s  G(s) 
s->-0 


K^,  Veloci-ty  Error  Constant.  The  velocity  error  constant,  is  defined  as 


K  =  lim  s  G(s) 
^  s-»-0 


Therefore, 


e(t)  =  R_ 

“  IV 
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vMch  is  the  error  in  displacement  (of  the  out^xit)  due  to  a  ramp  input. 
For  a  lype  0  ^stera 

K  =  lim  s  =  0, 

V  ft  0 

s-^0 


and  the  resulting  steai^-state  error  is  infinite  (Figure  13.43) 


c(t) 


00 


FIGURE  13.43.  STEADY-STATE  ERROR  -  TYPE 
"0"  SYSTEM,  RAMP  INPUT 
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FIGURE  13.44.  STEaDY-STftlE  RESPONSE  OF  A  TXPE  1 
SYSTEM  WITH  A  RAMP  INPUT 

For  2  and  3  systems  =  <»>  and  the  resulting  steai^-state 

error  is  zero. 

13.10.3  Parabolic  Input 

2 

Consider  the  input  r(t)  =  Rt 


e(t)  =  lim 
s^O 


=  lim 
s-»^0 


liin  s^  G(s) 
s-^0 


u_j(t)/2,  R(s)  =  R/s^ 


1  +  G(s) 

R 


s^  +  G(s) 
R 


K  ,  Acceleration  Error  Constant.  Hie  acceleration  error  constant,  K_,  is 

'  ■  ■  ■.■  '  -  S. 

defined  as 


k  =  lim  s^  G{s) 


Uie  steacty-state  error. 


e(t) 


ss 


is  the  error  in  displaceraent  (of  the  output)  due  to  an  acceleration  type 
ii^ait. 

For  a  lype  0  intern 


K  =  lim  K.  =  0 


and  for  a  Type  1  system 


s^  K, 


lim 

s-K) 


0 


For  lype  0  and  Type  1  ^tems  a  parabolic  input  vdll  result  in  a  pardsolic 
output  with  the  .  steacfy-state  error,  e(t)gg/  increasing  to  infinity  (Figure 
13.45) . 
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PARABOLIC 

INPUT 


C(t) 


FKURE  13.45.  STEADY-STATE  RESPCNSE  CF  TYPE  0  AND  TYPE 
SYSTEMS  TO  A  PARABOLIC  INPDT 


Etor  a  2  system 


s^  K, 


K_ 


Um 

s-»-0 


the  overall  gain  of  the  transfer  function. 

The  steac^-state  error,  is  ’equal  to  R/K.,  (a 

(Eigure  13.46) . 


constant) 


c(t) 


=  CONST. 


FIGURE  13.46.  STEADY-STATE  RESPCWSE  OF  A  TYPE  2 
SYSTEM  TO  A  PARABOLIC  INPOT 


For  a  I^pe  3  syston 


K  =  lim  — ^ 
®  s^O 


Therefore  the  stea^-state  error  is  zero. 

The  information  that  has  been  developed  is  presented  in  tabular  form  in 
Table  13.3. 
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TABLE  13.3 


Error  Constants 


System  I  Step  I  Banp  I  I^abolic  I  Step  Iiput 


Stea(^-State  Error 


Parabolic 

Bamp  Input  H^t 


=  K  =  \  =  ®<^)ss  =  -TTY^  ®<^Us  =  -r-  =  F 


13.10.4  Steady-State  Response  of  the  Control  Variables 

Ihe  foregoing  discussion  has  been  looking  at  the  steac^-state  error, 

based  on  a  specific  iiput  to  a  known  plant,  G(s) .  It  is  also  interesting  to 

look  at  tihe  stea^-state  value  of  the  control  variable,  c(t)  ,  for  a  known 

ss 

steac^-state  error  signal,  e(t) 

SS 

Consider  again  the  following  equation: 


Kn  (t^  s  +  1) +  1) 
(t^  s  +  1)  (Tj^  s  +  1) 


Rewriting  yields 


(t^s  +  1)  (Tj^S+  1)  ... 

^  (t^  s  +  1)  (t2  s  +  1) 


s^  C(s) 
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inlying  the  final  value  theorem 


e(t) 


ss 


lim  s  E(s) 
s-*>0 


lim 

s-^0 


*s  (t  s  +  1)  (Tj^  s  +  1)  . . . 
s  +  lMt2  S  +  1)  ... 


s"  C(s) 


lim  s 
s>0 


[s^  C(s)3 


Recall  the  differential  theorem 

[D"c(t)]  =  s"c(s) 


with  inificd  conditions  equal  to  zero. 

inlying  the  -Pinal  value  theorem  to  the  differential  -theorem  yields 

lim  s  [s*^  C(s)]  =  dP  c(t)gg 

s-^-0 


We  may  new  write 


d”  c(t) 


e(t) 


SS 


ss 


K 


n 


or 


From  this  equa-tion  and  the  characteristics  of  the  systems  as  diown  in 
Table  13.3,  the  following  conclusions  are  drawn  regarding  the  stea^-s-tate 
response; 


a.  A  type  0  sys-tem  is  one  in  vhich  a  cons-tant  actuating  signal 
a  constant  valuo  of  tho  output  ^ 
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b.  A  Typ©  1  syst63n  is  one  in  vAiich  a  constant  actuating  signal 
maintains  a  constant  rate  of  change  of  the  out^xit,  i.e. , 


c. 


e(t) 


ss 


D  c(t)33 


A  TVpe  2  systan  is  one  in  which  the  seocand  derivative  of  the  ^tput 
is  maintained  constant  by  a  constant  actuating  (error  signal)  i.e. , 


='«ss 


13.10.5  Deterniininq  System  Type  and  Gain  Fran  the  Bode  Plot 


System  type  and  gain  can  be  obtained  frcm  a  Bode  Plot  with  the  system  in 
unity  feedback  form.  The  slope  of  the  low  frequency  portion  of  the  In  curve 
determines  the  system  type:  a  0  db/dec  slope  represents  a  Type  0  system,  -20 
db/dec  a  Type  1  system,  and  a  -40  db/dec  slope  a  Type  2  system.  The  system 
gain  (K^)  is  determined  from  the  Im  plot  by  projecting  a  vertical  line  frcm 
0)  =  1.0  rad/sec  to  the  low  frequency  asymptote  (or  its  projection)  and  reading 
across  horizontally  the  Im  value.  This  Im  value  represents 

Im  G{ja))  =  20  log 


Depending  on  system  type,  =  K^,  or  K^.  Figure  13.47  illustrates  how 
the  Bode  Plot  is  used  in  the  manner  described. 
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FIGURE  13,47.  SYSTEM  TYPE  AT®  GAHT  FRC»1  A  BODE  PLOT 


13.87 


13.10.6  Suttnary 

TSie  static  error  constant  can  be  used  to  quickly  determine  the  ability  of 
control  system  to  follow  a  specific  input.  The  application  of  error  constants 
is  not  1  to  ^tems  with  inputs  classifi^  as  one  of  the  three  basic 
types  of  test  signals.  Ibr  linear  systems,  the  concept  can  easily  be  extended 
to  systems  with  inputs  that  can  be  represented  by  a  polyncmial,  i.e.. 


r(t) 


TSie  stea(^-state  error  is 


+  t  + 


(t) 


e(t) 


ss 


R 


a  superposition  of  the  errors  due  to  each  input  signal  ccrponent  acting  alone. 

The  chief  advantage  to  the  foregoing  approach  to  steady-state  response  is 
the  ease  and  timeliness  of  arriving  at  the  ansvrer.  The  chief  disadvantage  of 
tile  error  constant  approach  is  that  only  one  of  the  constants  has  a  finite 
value  vjhich  is  not  zero  or  infinity  for  a  particular  n-type  system.  In  cases 
vdiere  the  stea^^-state  error  is  a  function  of  time,  the  error  constant 
approach  only  gives  an  answer  of  infinity  for  e(t)gg  and  does  not  provide  an 
indication  of  how  the  error  varies  with  time.  Even  though  the  steady-state 
error  may  turn  out  to  be  infinite,  for  an  actual  problem,  the  injxit  may  be 
applied  for  a  finite  time,  thus  the  error  will  be  finite.  This  finite  error 
may  be  well  within  the  specifications. 

It  iiould  appear  desirable  to  select  a  large  value  of  the  overall  gain 
of  the  transfer  function,  to  minimize  the  steady-state  error;  but  not  without 
a  penalty.  Itoo  large  a  value  of  K  may  force  the  ^stem  unstable.  As  \^/e  will 
see  when  we  get  to  root  locus  analysis,  an  adjustment  of  the  system  gain 
effects  both  the  natural  frequency  and  the  danping  ratio  for  a  closed-locp 
system  with  ccnplex  poles.  In  many  cases  the  exact  value  of  vdiich  results 
in  unstable  system  operation  may  be  found  by  analysis.  Routh's  criterion,  for 
exanple,  will  provide  the  value  and  the  application  of  the  criterion  may  be 
fejund  in  the  literature  (Reference  1). 
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13.11  ROOT  LOCUS 


An  accurate  prediction  of  the  system^  s  performance  can  be  obtained  by 
deriving  the  differential  equation  of  a  control  system  and  then  determining 
its  solution.  This  approach  is  not  feasible,  however,  for  any  but  the 
simplest  system.  Not  only  is  the  direct  solution  method  extremely  tedious, 
but  If  the  response  does  not  meet  the  required  specifications,  no  indication 
is  given  of  hov  to  ivaprove  its  performance. 

The  aim  of  the  design  engineer  is  to  predict  the  performance  of  the 
system  without  solving  its  equations  of  motion.  Also,  he  would  like  the 
analysis  to  indicate  how  to  modify  the  system  in  order  to  produce  the  desired 
response  characteristics.  Several  methods  are  available  which  both  predict 
stabili^  and  indicate  the  of  condensation  required.  Of  those,  the  Bode 
plot  has  alreacty  been  discussed  and  the  root  locus  will  be  discussed  next. 
Another  technique  is  Nyquist  criterion.  The  theory  and  application  of  root 
locus  will  be  described. 

Definition;  Olie  root  locus  is  a  plot  of  the  roots  of  the  characteristic 
equation  of  the  closed-loop  system  as  the  gain  is  varied  from  zero  to 
infinity.  The  definition  itself  presents  the  underlying  theory  of  the  root 
locus  method.  The  primary  objective  is  to  determine  system  stability.  This 
leads  to  another  question.  What  determines  stability?  The  answer,  is  the 
transient  solution,  vhich  is  determined  from  the  roots  of  the  characteristic 
equation,  vhich  cannot  have  positive  real  parts  and  be  stable. 

The  general  approach  used  in  the  development  of  the  root  locus  technique 
will  be  to  plot  a  root  locus  for  a  sinple  system  the  hard  way,  i.e., 
successive  analytic  solutions  for  the  roots  of  the  characteristic  equation  for 
selected  values  of  static  loop  gain,  K.  Then  the  significance  of  the  root 
locus  will  be  discussed  for  the  sinple  system,  and  then  for  any  system. 
Lastly,  some  rules  will  be  developed  which  permit  quick  plots  to  be  drawn 
using  relatively  little  labor. 

13.11.1  Poles  and  Zeros 

This  section  will  define  vhat  is  meant  by  poles  and  zeros  and  also 
discuss  their  relationship  in  functions  of  interest  here.  Consider  the  system 
represented  by  the  block  diagram  in  Figure  13.48. 
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FIGURE  13.48.  CDOSED-LOOP  SYSTEM 


\4iere 

C(s)  ^  G(s) 

R(s)  1  +  GH(s) 

and  viiere  in  general 

G(s) 

and 

H(s) 

viiere  the  numbers  •••;  •••'  ***'  "* 

real,  oonplex  or  zero. 

We  now  define  tsro  nev  terms; 

ZERO  —  A  zero  of  a  function  (like  G(s) )  is  a  value  of  s  that  makes  that 
function  zero. 

POTF.  —  A  pole  of  a  function  (like  G(s) )  is  a  value  of  s  that  makes  that 
fimction  go  to  infinity. 


K  (t.s  +  1) (t-s  +  1) (...) 
n  1  ^  _  _ 


n  (t  s  +  1)  (t,  s  +  !)(...) 
S  U  ^ 


=  Kj^G  (s)  = 


Dg(s) 


(13.71) 


(t  s  +  1)  (t  s  +  1)  (...) 

V  Q  g _ 

n  (t-S  +  1) (TjjS  +  1) (...) 


=  (s)  = 


K 


(13.72) 
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Bar  Exaitple: 


s  =  -  —  is  a  zero  of  G(s) 


s  =  -  ~  is  a  zero  of  H(s) 

T 

a 


s  =  -  —  is  a  pole  of  G(s) 

T 

a 


s  =  -  —  is  a  pole  of  H(s) 

“'ll 

In  terms  of  the  s-plane,  (Figure  13.49) ,  this  means  that  there  are  values 
that  cause  G(s)  and  H(s)  and  incidentally  their  product  G(s)  H(s)  to  be  zero 
and  to  be  infinite.  Figure  13.50  is  a  plot  of  the  function  G(s)  H(s) . 


FIGURE  13.49.  s-PIANE  FIGURE  13.50.  SURfftCE  CF  G(s)  H(s) 
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The  value  of  s  which  results  in  an  infinite  value  of  G(s)H(s)  is  the  pole 
of  G(s)H(s). 

The  pole  gets  its  name  frcm  the  appearance  that  a  graph  of  the  magnitude 
of  G(s)H(s)  makes  as  "s"  assumes  values  near  the  point  +  13.51). 


FIGURE  13.51.  A  POiLE  OF  G(s)H(s) 


Note  that  the  poles  and  zeros  of  the  function  G(s)H(s)  corpletely 
describe  the  fimction.  When  we  take  the  conposite  function  like 
G(s)/[1  +  G(s)H(s)],  vhere  presumably  we  know  the  poles  and  zeros  of  G(s)  and 
H(s),  one  must  effircise  caution  regarding  the  transference  of  this  information 
to  the  ccnposite  function. 


Ibr  Example: 

A  zero  of 
A  pole  of 
A  zero  of 
A  pole  of 
A  zero  of 


G(s)  is  also  a  zero  of  G(s)/[1  +  G(s)H(s)] 

G(s)  does  not  resxilt  in  a  pole  or  zero  of  G(s)/[1 
H(s)  does  not  result  in  a  pole  or  zero  of  G{s)/[1 
H(s)  is  a  zero  of  G(s)/Il  +  G(s)H(s)] 

1  +  G(s)H(s)  is  a  pole  of  G(s)/[1  +  G(s)H(s)] 


+  G(s)H(s)] 
+  G(s)H(s)] 
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NOW,  since  we  want  to  see  if  the  transients  die  out  let  us  take.,  the 
expression  C{s)/R(s)  (Equation  13.69),  solve  it  for  C(s)  and  assume  some  form 
of  excitation  R(s).  Actxially  ai^  form  of  excitation  (sine,  unit  step,  tmit 
vansp,  etc. )  may  be  used. 

Substituting  Equations  13.71  and  13.72  into  Equation  13.69  yields 

C(s)  „  G(s) 

ITiT  1  +  G(s)H(s) 


K 


'G 

TpsT 


Tpi) 


\  N^(S)  D„(S) 


D^(S)  D^(S)  +  N^(S)I^(S) 


(13.73) 


The  zeros  of  Dq(s)  Dg(s)  +  Nq(s)  Ng(s)  are  the  same  as  the  zeros  of 
1  +  G(s)H(s)  and  Equation  13.73  can  be  factored  into  the  form 


Kn  Ne(5)  D,(s) 


K„  N,(S)  D,(S) 


C(S) 


(s  -  rj^  )(s  -  rj  )(s  -  r^ )  ...  (roots  from  Root  Locus) 


(13.74) 


>4iere,  for  convenience  we  let  0^  (s)  =  1,  the  unit  inpilse  function.  By 
partial  fractions  Equation  13.74  can  be  expanded  into  the  form 


C(s)  = 


s  -  r. 


s  -  r. 


s  -  r. 


vAiere  r^  are  the  zeros  of  1  +  G(s)H(s)  and  the  poles  of  G(s)/(1  +  G(s)H(s)]. 
The  factors,  r^ ,  may  be  real  or  conplex  and  positive  or  negative.  Note  that 
the  inverse  transform  of  each  element  leads  to  an  exponential  term.  Assuming 
r^  is  real  and  positive,  then  (r^  =  +  )  and 
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s 


s 


(13.75 


+a,  t 
1 

e 


UNSTABLE 


and  if  is  negative  then  ( ) 


As  -»st 


STABLE 


(13.76), 


In  the  first  case  (Equation  13.75),  the  anplitude  of  the  transient  term 
gets  large  as  time  gets  large  because  of  the  term. 

In  the  second  case  (Equation  13.76),  the  transient  term  disappears 

-a.t 

because  as  time  gets  large  e  goes  to  zero. 

Thus,  if  a  system  under  investigation  has  any  positive  real  poles  of 
G(s)/[1  -f  G(s)H(s)]  or  a  positive  real  zero  of  1  +  G(s)H(s)  then  the  system  is 
unstable. 

Conversely,  if  the  system  being  investigated  has  all  negative  real  poles 
of  G(s)/(1  +  G(s)H(s)]  or  all  negative  real  zeros  of  1  +  G(s)H(s),  then  the* 
system  is  stable. 

If  we  assume  that  r^  is  conplex  then  we  know  there  exists  another  zero  of 
1  +  G(s)H(s)  \diich  is  the  conplex  conjugate  of  r^,  namely . 

Ttiis  pair  of  zeros  of  1  +  G(s)H(s)  leads  to  a  term  in  the  partial 
fraction  expansion  \diere  r^  *  and  -  jw^  of  the  form 

(s  -  oj"  +  (wj" 

\diich  has  the  inverse  transform  of  the  form 

K  e“c^  cos  (w^t  +  4)^) 

\diere  if  is  positive,  =  +(t^  ,  then  we  get  an  exponentially  increasing 
cosine  term  (Figure  13.52). 
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-  K  •  +  0c) 


,:"c^ 


f  C*  FORMS  THE  ENVELOPE 


FIGORE  13.52.  EXPOilENTIALLy  INC3?EASING  C30SINE  TERM 


HCxrever,  if  a  is  negative,  a  =  -o  .  Then  the  response  is  of  the  form 


Ife  oos  (a)^t  +  4'^)  \Aiich  is  a  cosine  term  with  an  envelope  that  decreases 
with  time  (Figure  13.53) . 


Ke  "C‘co8(ajj.t  +  6_) 


isra 


FORMS  THE  ENVELOPE 


FIGDRE  13.53.  EXPC»IINTIALLy  DBC3?EASING  COSINE  TEFM 

Thus,  we  conclude  that  if  a  ocraplex  sero  of  1  +  G{s)H(s)  has  a  positive 

real  part,  a  =  +a^,  then  the  system  is  unstable  and  if  a  cotplex  zero  of 

1  +  G(s)  H(s)  has  a  negative  real  part,  a  -  -a  ,  then  the  system  is  stable. 

c 

Actually  the  conditions  for  ireal  zeros  and  cotplex  zeros  are  the  same: 
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REAL  PART  POSITIVE  —  SYSTEM  UNSTABLE 
REAL  PART  NEGATIVE  —  SYSTEM  STABLE. 

Now  v^t  is  the  significance  of  the  location  of  the  zeros  of  1  +  G(s)H(s) 
x^wn  the  s-plane?  Looking  at  the  s-plane  we  find  that  if  ANY  zeros  of 
1  +  G(s)H(s)  are  in  the  RHP  the  system  is  xinstable.  If  ALL  zeros  of 
1  ^  G(6)H(8)  are  in  the  LHP  the  system  is  stable. 

Now  knowing  that  instability  is  caused  by  a  zero  or  zeros  of 
1  +  G(s)H(6)  with  a  positive  real  part,  the  problem  of  determining  stability 
degenerates  to  the  problem  of  determining  vdiether  or  not  there  are  any  zeros 
of  1  +  G(s)H(s)  in  the  RHP  or,  equivalently,  \diether  1  +  G(s)H(s)  does  indeed 
have  a  zero  or  zeros  with  positive  real  parts. 

13.11.2  Direct  Locus  Plotting 

The  example  to  be  used  for  direct  root  locus  plotting  will  be  the  second- 
order  Intern  tdiose  differential  equation  and  transfer  functions  were  derived 
earlier.  The  system  is  shown  in  Figure  13.54  and  Equation  13.67  gives  the 
forward  transfer  function:  H(s)  <=  1. 


FIGURE  13.54.  UNIT  FEEDBACK  SYSTEM 
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G(s) 


P 


Is^  +  bs 

The  following  values  will  be  assumed  for  the  constants 

1  =  1 
b  =  2 
V  =  K 


Equation  13.67  becones 


G(s) 


K 

s(s  +  2) 


(13.77) 


The  problem  is  to  determine  the  roots  of  the  charac±eristic  equation  for 
all  values  of  K  and  to  plot  these  roots  in  the  s-plane. 

Fran  Equation  13.69  the  system  closed-loop  transfer  function  is 


G(s) 

1  +  (3i{s) 


K 

s(s  2)  ^  K 


Ihe  system  characteristic  equation  is 

s^  +  2s  +  K  =  0 


(13.78) 


The  roots  of  Equation  13.78  are 


s^^2  =  -  1  ±  Vl  -  K 

The  location  of  roots  for  various  values  of  K  is  shown  in  Table  13.4. 
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TABLE  13.4 

CLDSEEHLOOP  BOOT  LOCATIONS  AS  A  FUNCTION  OF  K 


K 

®1 

0 

0  +  jO 

“2  -  jO  (open-loop  poles) 

1/2 

-.3  +  jO 

-1.7  -  jO 
/ 

1 

-  1  +  jO 

-1  -  jO 

2 

-  1  +  jl 

-1  -  jl 

3 

-  1  +  j  VT 

-1  -  j  VF 

00 

-  1  +  j® 

-1  -  j® 

The  points  fron  Table  13.4  are  plotted  in  the  s-plane  of  Figure  13.55 
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The  root  locus  of  Figure  13.55,  ^ich  is  the  locus  of  roots  of*,  the 
characteristic  equation  as  a  function  of  gain,  quickly  indicates  vdiether  or^ 
not  the  system  is  stable,  and,  also,  the  form  of  the  transient  response  for 
any  selected  value  of  K.  From  the  plot,  it  can  be  seen  that  for  0  <  K  <  1  the 
roots  are  real  and  negative  resulting  in  ex^xmential  decay  from  each  root. 
For  1  <  K  <  «»,  however,  the  roots  are  conplex  with  the  real  part  negative. 
The  corresponding  transient  response  is  oscillatory  within  an  exponentially 
decaying  envelope.  For  exanple,  for  K  1.5  the  approximate  values  of  s  from 
the  locus  are 


Sj  =  -1  +  j  0.5 
Sj  «  -1  -  j  0.5 

The  system  transient  response  for  K  «  1.5  is  then 

c(t)  =  Ce~^  cos(0.5t  +  ♦) 

The  time  involved  in  constructing  a  root  locus  for  a  conplex  system  in  this 
manner  is  obviously  prohibitive.  This  difficulty  will  be  overcome  later.  | 
To  further  discuss  the  significance  of  the  s-plane  and  the  root  locus  we 
will  consider  a  second-order  characteristic  equation  in  its  standard  form 

s^  >f  2C<«)_s  +  »_^  «  0 

n  n 


The  roots  of  this  equation  are 


ff  (0^ 


(13.79) 


In  order  to  realize  vhat  this  means  in  the  s-plane,  refer  to  Figure  13.56a. 
Any  arbitrary  value  of  s  will  have,  from  Equation  13.79,  a  real  part  a  =  Cw^ 
and  an  imaginary  part  ^  1  -  .  From  these  values  shown  in  Figure 
13.56,  and  the  Pythagorean  theorem,  the  magnitude  of  the  position  vector,  s  is 
found  to  be  equal  co^ .  The  angle  <j>  is  also  significant  since 
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CX>S  (|) 


Figure  13.56a  summarizes  this  information  and  shews  hew  parameters  important 
to  the  transient  response  can  be  easily  obtained  from  the  position  of  roots  in 
the  s-plane.  From  the  root  locus,  then,  the  transient  response 
characteristics  for  all  values  of  gain,  K,  can  be  seen  at  a  glance. 


13.11.3  Angle  and  Magnitude  Conditions 

New  that  the  value  of  the  root  locus  has  been  established,  the  rules  will 
be  developed  vhich  permit  simplified  plotting  of  ccn^lex  ^tans.  Ihese  rules 
are  based  on  two  conditions,  the  angle  condition  and  the  magnitude  condition 
vhich  evolve  &om  the  characteristic  equation. 

As  stated  many  times,  the  closed-loop  control  system  may  be  represented 
by 


G(s) 

1  +  GH(s) 


from  vhich  the  characteristio  equation  is 


1  +  GH{s)  =  0 
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Bie  ^tero  open-loop  transfer  function  can  be  writtai  in  the  following  form 


GH(s) 


K(s  -  z^)  {s  ~  z^)  ,  ,  . 
(s  -  pj)  (s  -  P2)  .  .  . 


(13.80) 


\diere  the  static  gain,  K,  is  factored  out  and  the  z's  and  p's  are  the 
open-loc^  zeros  and  poles  respectively. 

Since  the  values  of  s  to  be  determined  must  satisfy  the  relationship 


1  +  GH(s)  =  0 


we  can  say 


GH(s)  =  -1 


(13.81) 


but,  since  s  is,  in  general,  ccnplex,  GH(s)  is  then  a  ftmction  of  a  coiplex 
variable  and  Equation  13.81  can  be  written 


GH(s)  =  iej(l  +  2n)ii  ^  ^  /(l  -i-  2n)Tr 


(13.82) 


where 

n  =  0,  +  1,  +  2,  +  3,  .  .  .  . 

Equation  13.82  says  that  in  order  for  the  value  of  s  to  be  a  z^o  of  1  +  GH(s) 
the  magnitude  of  the  coiplex  quantity  GH{s)  must  be  equal  to  1  and  the 
argument  be  sane  odd  multiple  of  ir.  Hence  the 

Magnitude  condition 

|(3i(s)|  =  1  (13.83) 

Angle  condition 

/gH(s)  =  (l  +  2n)Tr  (13.84) 

n  =  0,  +  1,  +  2,  +  3  .  .  . 
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Substituting  Equation  13.80  into  Equation  13.84,  the  angle  condition  beoanes 


lC{s  *  ”  ^2^  *  *  * 

(s  -  p, )  (s  -  p,)  .  .  . 


=  (1  +  2n)  ir 


(13.85) 


pjwffta  each  fecbor  of  C3i{s)  can  be  represented  by  a  vector  in  the  silane  frcm 
the  pole  or  zero  to  the  s-point  in  question,  Equation  13.85  can  be  written 

A  .  (1  .  2n), 

/s  -  pl-^  A  -  P2  • 

or 

/s  -  +  /s  -  +  ..."  /s  -  P^  -  A  ~  ^2  ~  *  •  *  (1  +  2n)Ti  (13.86) 

Thus,  »ging  the  angle  condition,  any  point  in  the  s-plane  can  be 
investigated  to  determine  vhether  or  not  it  is  a  point  on  tiie  root  locus  by 
measuring  the  angles  of  the  vectors  from  the  poles  and  zeros  to  the  point  in 
question,  and  adding  them  according  to  the  left  side  of  Equation  13.86.  If 

this  sum  equals  an  odd  multiple  of  ir,  the  value  of  s  satisfies  the 

characteristic  equation,  and  is  on  the  root  locus. 

Figure  13.57  demonstrates  the  application  of  the  angle  condition.  The 
vectors  representing  GH(s)  for  s  =  s^  are  shown  in  Figure  13.57.  The  angle 

condition  test  to  determine  if  Sq  is  on  the  root  locus  is 

~  ®1  “  ®2  ”  (1  +  2n)Tr  (13.87) 


From  the  figure. 


^  =  170° 

=  195° 
02  =  140° 
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OherefDre,  Sq  is  not  on  the  locus  but  is  fairly  close.  Successive  tries  will 
allow  converging  on  the  point  that  satisfies  the  characteristic  equation. 
CJotibining  Equation  13.80  and  13.83  the  magnitude  condition  becanes 


K  •  S  ■“  2^^ 

• 

s  - 

•^2  • 

•  • 

r' 

'ifi 

s  - 

s  - 

•  P21  •  . 

•  1=  -  P21 

•  •  • 

(13.88) 

K  =  - 

s  - 

! 

• 

s  - 

1  •  •  •  • 

Equation  13.88  says  that  the  magnitude  of  all  the  vectors  from  the  poles  to 
the  s  point  in  question,  divided  by  the  magnitudes  of  the  vectors  frcm  the 
zeroes  is  equal  to  the  gain,  K.  This  condition  allows  us  to  determine  the 
value  of  gain  frcm  ary  point  on  a  loais  in  the  s-plane. 

•She  reader  is  encouraged  to  check  the  example  system's  root  locus 
(Figure  13.55) ,  to  ensure  it  is  consistent  with  both  the  angle  and  magnitude 
condition. 


13.11.4  Rules  for  Boot  Locus  Construction 

The  following  rules  f or  K  >  0  will  allow  a  sketch  of  the  root  locus  to  be 
drawn  quickly.  Ihese  rules  are  based  upon  the  angle  condition  and  an  analysis 
of  the  characteristic  equation. 

a.  Ihe  nuTiber  of  branches  of  the  locus  is  equal  to  the  number  of 
open— loop  poles  (i.e. ,  the  order  of  the  characteristic  equation) . 


1  +  GH(s) 


K(s  -  z^)  ...  (s  “  z^) 
(s  “  ^1^  ...  (s  “ 


(13.89) 


(s  -  Pj)  .  .  .  (s  -  Pp)  +  K  (s  -  z^)  .  .  .  (s  -  Zg)  -  0 


Since  we  have  assumed  a  rational  polynomial,  P  ^  Z,  and  the  hipest 
order  of  s  is  P. 

P  is  the  of  open  loop  poles,  Z  is  the  number  of  opoi  loop 

zeros. 


b.  Tlie  loci  branches  b^in  at  the  open-lobp  poles  viiere  K  =  0.  If  ws 
write  the  characteristic  equation  with  the  static  gain  factored  out 


1  +  GH(s)  =  1  +  KGH'  (s)  =  0 

then 


KGH'  (s)  =  -  1 


GH'  (s)  =  -  1/K  (13.90) 


for  K  =  0,  GH'  (s)  =  <*>,  vtoch  means  s  is  at  a  pole  of  GH(s) . 

c,  Ihe  branches  end  at  the  open-loop  zeros  viiere  K  =  «.  Eron  Equation 
13.90,  vAien  K  =  »,  GH'  (s)  =  0,  viiich  is  a  zero  of  GH(s) .  When  P  >  Z 
(the  usual  case)  the  additional  branches  end  at  s  =  <»  vdiich  may  also 
be  considered  an  open-loop  zero. 

d.  Bie  loci  branches  that  do  approach  s  =  «*>  do  so  asymptotically  to 
radial  lines  centered  at 


_  E  Rs  Cp's3  -  E  He  fz's  3 
P  -  Z 


e.  The  angles  of  these  asynoptotes  are  given  by 

_  (1  +  2n)Tr 

'  P  -  Z 


vhere 


n  =  0,  +1,  +  2... 


f .  A  point  on  the  real  axis  is  on  a  locus  branch  if  it  is  to  the  left 

of  an  odd  number  of  open-loop  poles  and  zeros.  This  can  be  seen 
from  Figure  13.58  in  vhich  it  is  obvious  that  the  net  angular 
contribution  from  a  pair  of  ccnplex  conjugate  poles  or  zeros  to  a 
search  point  on  the  real  axis  is  0^.  The  real  axis  poles  and  zeros 
contribute  zero  or  ir  vhen  the  search  point  is  to  the  ri^t  or  left 
of  the  pole  or  zero  respectively.  Therefore,  the  angle  condition  is 
satisfied  only  to  the  left  of  an  odd  number  of  open-loop  poles  and 
zeros. 
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FISURE  13.58.  RE^L  AXIS  LOCI 


g.  Itie  angle  of  departure  approach  of  a  loci  branch  from  (to)  a  ccnplex 
pole  (zero)  can  be  found  vising  the  angle  condition  as  follows.  If  a 
search  point  is  chosen  e  distance  from  a  coiplex  pole  (zero) ,  the 
angle  of  this  vector  is  the  departure  (approach)  angle  at  that  pole 
(zero) . 

Consider  the  situation  of  Figure  13.59  vdiere  it  is  desired  to  find  the 
departvire  angle  from  p^.  If  the  magnitude  of  the  vector  fron  p^  is  e  then  the 
angle  of  the  vectors  from  the  other  poles  and  zeros  can  be  measured  directly 
to  p^.  Solving  the  angle  condition,  Bguation  13.86  for  ~  (an  vinkncwn 
quantity)  will  yield  the  departvire  angle. 
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FIGUKE  13.59.  DEPARIUKE  ANGLE  DETERMINATION 
Ecjuation  13.86  says 

y3-.-/s-p^-/s-p,-/s-P3-  /s-p,  ^  U.2n), 

Ercn  Figure  13.59 

((•  ”  6^  “  ©2  ”  ^2  ~  (1  2n)  TT 

(j)  =  45° 

02  =  90° 
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155® 


*  45®  -  90®  -  135®  -  155®  -  n 

-  -  515®  »  -  155® 

h.  The  points  at  iidiich  the  loci  branches  leave  or  enter  the  real  axis 
are  sometiines  called  "breakaway"  or  "breakin"  points /  respectively. 

Figure  13.60  illustrates  the  conpitation  of  a  breakaway  point  for  the 
system  v^ose  open-loop  transfer  function  is 

*  s{s  +  IHs  +  2) 

From  the  characteristic  equation, 

_ K  , 

s(s  +  l)(s  +  2) 

K  «  -  s  (s  +  l)(s  +  2) 

.  K  *  -  s^  -  3s^  -  2s 

To  firrf  the  point  vdiere  K  is  meucimum,  we  can  differentiate  the  expression  for 
K  and  set  it  equal  to  zero.  Ihe  solution  of  this  equation  should  produce  the 
desired  restilt. 


e. 


e. 


Departure  angle  is 


e. 


g  =  -  3s^  -  6s  -  2  =  0 

Sj  2  =  -  1  +  0.574 

s  =  -  1.574  or  -  0.426 
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4 


S-PLANE 


FIGUPE  13.60.  BREAKmY  POINT  COMPUEATION 
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It  is  obvious  from  Figure  13.60  that  the  latter  solution,  s  *  -0.4257.,  is 
the  meaningful  answer  beca\ise  it  ai^ars  on  a  real  axis  locus.  This,  then  id| 
the  breakaway  point. 

This  method  of  computing  breakaway  points  is  restrictive  in  that  more 
conplex  systems  with  higher-order  transfer  functions  are  difficult  to  solve 
since  dVds  is  higher  than  second-order. 

The  fact  that  the  root  locus  is  always  symmetrical  about  the  real  axis  is 
advantageous  in  that  only  the  tqc^r  half  of  the  s-plane  need  be  plotted.  The 
lower  half  is  just  a  mirror  image.  The  reason  for  this  should  be  obvious  xdien 
one  considers  that  ccmplex  roots  always  appear  in  conjugate  pairs. 

The  rules  above  will  permit  a  sketch  of  the  root  locus  very  quickly.  If 
more  accurate  data  is  required,  the  branches  can  be  checked  using  the  angle 
conditicm  and  a  protractor.  After  the  locus  is  conplete,  the  values  of  K  that 
are  deemed  inportant  can  be  conpited  using  the  magnitude  condition. 

Boot  Locus  Exanples: 

Exanple  1: 


“  s(6.5s  +  l)(0.is  +  1) 

H(s)  -  1 

Rewriting  in  the  more  useful  form 

'  s(s  +"2?fr+?] 

Letting  lOK'  »  K,  we  can  plot  the  locus  and  calibrate  it  in  gain.  In  the 
actual  system,  however,  the  gain  selected,  K',  will  be  a  factor  of  10  less  than 
that  found  from  the  locus  plot. 

From  the  previous  equation 


Pi 


0,  Pj  =  -2,  P3  =  -5 


P  =3,  Z  =  0 
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J^plying  tiie  rules  developed  in  the  previous  section/  the  root  locus  of  Figure 
13.61  is  plotted. 

a.  Bie  number  of  branches  is  3. 

b.  The  locus  branches  begin  at  the  poles  s  =  0#  -2,  -5,  vAiere  K  =  0. 

c.  Since  there  are  no  open-loop  zeros,  the  3  branches  will  end  along 
asymptotes  \diose  real  axis  intercepts  are 


a 


E  Be  [p's]  -  Z  Be  [z's3 
P  -  Z 


a 


[0  -  2  -  51  -  [0] 
3T-0 


a  =  -  7/3 

The  asynptotes  angles  are 

_  (1  +  2n)Tr 
7  P  -  Z 

y  =  (1  +^2n)TT  ^  gpO^  j^gQO^  3qqO 

d.  Real  axis  loci  exist  between  s  =  0  and  s  =  -2,  and  frcm  s  =  -5  to  s 
=  “<*>. 

e.  There  are  no  ccmplex  poles  or  zeros  so  rule  (g)  does  not  apply. 

f.  The  breakaway  point  is  found  as  follows: 
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K-70 

- h- 
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GH(S} 


-1 


K 

s(s  +  2)(s  +  5) 

-  s(s  +  2)(s  +  5) 

-  -  7s^  -  10s 

-  3s^  -  14s  -  10  =  0 

-  3.79,  -  0.88 

Since  -  3.79  is  not  on  a  locus  branch,  s  «  -0.88  is  the  breakaway  point. 

Now  that  the  root  locus  is  plotted,  the  desired  values  of  K  can  be  found. 
It  is  generally  tiseful  to  know  the  value  of  K  for  \duch  the  system  becomes 
unstable.  By  measuring  the  length  of  the  vectors  from  each  of  the  poles  to 
the  point  v4iere  the  locus  crosses  the  imaginary  axis  into  the  LHP,  we  find 

K  -  (3.16)(3.7)(6.0)  -  70 

Bie  system  is  neutrally  stable  at  K  =  70  and  the  frequency  of  the  oscillation 
will  be  3.16  ra^sec.  For  K  >  70,  the  roots  of  the  characteristic  equation 
are  in  the  RHP  vdiere  they  have  positive  real  parts  and  the  system  is  unstable. 

If  some  specific  transient  response  parameter  is  required,  such  as 
C  “  0.5,  the  corresponding  vcilue  of  K  can  be  determined. 

Since 


K  » 
K  « 
dK 

ai ' 

s  ■= 


cos“^  t  ^ 
cos'^0.5  =  60“ 

Where  the  radial  from  the  origin  60“  from  the  negative  real  axis  crosses  the 
root  locus,  the  system  is  at  the  required  operating  condition,  and  the  gain 
for  this  point  is 
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K 


(4.4)  (1.7) (1.4) 


10.6 


To  find  the  total  transient  response,  however,  it  is  necessary  to  find  the 
point  on  all  branches  where  K  =  10.6.  Each  branch  contributes  one  to  the 
transient  re^nse  and  th^  all  must  have  the  same  value  of  K.  Because  the 
lower  half  plane  is  a  mirror  image,  the  lower  root  is  the  oonplex  conjugate  of 
the  i^per  value.  Bie  point  on  the  third  branch  along  the  real  axis  vdiere  K  = 
10.6  is  found  throu^  trial  and  error.  The  transient  solution,  then,  for  K  = 
10.6,  ?  =  0.5,  has  a  quickly  danped  pure  e}qx)nential  term  and  a  dominant  more 
slcwly  danped  ejqonentially  decaying  oscillation  at  =  1.4. 


®0<« transient  =  ««  '1-2^  +  ♦> 


Example  2: 


GH(s)  = 


K(s  +  2) 


s(s  +  3)  (s*^  +  2s  +  2) 


GH(s)  = 


K(s  +  2) 


s(s  +  3) (s  +  1  +  j) (s  +  1  -  j) 


inlying  the  rules: 


a. 

The  number  of  branches  is  4. 

b. 

The  locus  brandies  begin  on  the  open-loop 

"1  ±  j  • 

poles  at  s  =  0,  -3, 

c. 

Oie  branch  will  end  on  the  open  locp  zero  at 
3  branches  will  end  at  s  =  “  along  asymptotes 

s  =  -2.  The  remaining 
centered  at 

_  _  [0  -  3  -  1  -  1]  -  [  -  21 

4-1 

a  =  -  1 
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wi-di  asynptote  angles 


V  =  (1  +  2n)Tr 

'  4-1 

Y  =  60°,  180°,  300° 

d.  Iteed  axis  loci  exist  betwcjen  s  =  0  and  -2,  and  form  s  =  -3  to  - 

e.  Ihe  departure  angle  from  the  ccttplex  poles  is  found  by 

0  =  45°  -  135°  -  90°  -26.5°  +  180° 
e  =  -  26.5° 

Once  the  value  of  s  for  the  imaginary  axis  crossing  is  foiand,  a  graphical 
solution  for  K  using  the  magnitude  condition  can  be  used,  i.e. , 

K  =  (3.41)  (2.81)  (1.60)  (1.19)  ^  7  q. 


We  will  continue  this  development  to  determine  the  closed-loop  response 
for  ?  =  .5.  As  the  open-loop  transfer  function  was  specified 


GH 


_ K  (s  +  2) 

s(s  +  3)  {s^  +  2s  +  2) 


we  shall  further  specify  that 


G(S) 

H(S) 


K 


s(s^  +  2s  +  2) 


(s  +  2) 
(s  +  3) 


and 


It  can  be  seen  fron  the  root  locus  plot  that  the  transient  response  will 
have  two  pure  exponentially  decaying  terms  from  the  real  axis  branches  and  one 
oscillatory  damped  term  frcra  the  ccnplex  loci.  It  shoulci  also  be  noted  that 
the  maximum  c  possible  for  the  oscillatory  terms. is  c  =  .7  established  by  the 
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ppei>-lopp  oonplex  poles  where  K  =  0.  Althou^  the  C  =  .5  radial  crosses  the 
carplex  locus  in  the  vicinity  of  s^^  2  “  0.5  +  .'Sj,  the  exact  crossing  cein  be 
found  using  the  angle  condition.  A  point  on  the  c  =  >5  radial  near  the 
sketclted  ?.<  c\:s  is  tried  as  a  search  point  and  if  the  point  satisfies  the  angle 
condition,  the  point  is  on  the  locus  and  pins  down  the  exact  roots  in 
question.  A  couple  of  search  point  narrows  our  locus  to  s  =  -.56  +  .96j.  At 
this  point  the  angle  condition  is  satisfied,  i.e. ,  ys  -  -  p^ 

-  /s  -  p,  -  /s  -  p^  -  -  Pt  =  ir,  or  33.6°  -  120°  -  7^  -  21.2°  -  (-5^) 

=  179.6^  £  u. 

Next,  the  magnitude  condition  used  at  this  point  is  =  -.56  +  j.96,  to 
determine  the  value  of  gain,  K. 

K  =  is|  Is  +  3|  is+  1  til  Is  +  1  -ii 

|s  +  2| 

_  (2.36)  (2.025)  (1.12)  (.45) 

1.74 


=  1.532. 


So  far  we  have  determined  two  roots  (a  ccnplex  pair)  and  the  value  of  gain  for 
system  operation.  7V?d  other  roots  lie  on  the  real  axis  loci  (one  on  each 
branch)  and  are  found  by  trial  and  error  by  using  the  magnitude  condition  and 
finding  a  search  point  that  will  yield  a  value  of  K  =  1.532. 

A  search  point  at  s  =  -.8  is  determined  to  be  a  root  corresponding  to  a 
value  of  K  =  1.532,  i.e. , 


(.8)  (2.2) (1.025)^  _  T  c.  _  V 

(1.2)  “  ~  ^ 

(close  enoui^  for  a  graphical  soluticai) 

toother  search  point  at  s  =  -3.81  also  satisfies  the  magnitude  condition, 

i.e. , 


(.81)  (3.81)  (2.99)^ 
1.81 


1.525  =  K 


(close  enou^  for  a  graphical  solution) 
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13.12  CCWPENSATION  TECHNIQUES 

You  have  been  introduced  to  theory  and  technique  used  to  set  up  a  control 
system  problem  and  to  analyze  the  results.  Unfortunately  the  system  often 
needs  to  be  "fixed"  to  meet  the  desired  performance  specifications.  Biis 
"fixing"  is  called  conpensation  and  is  used  to  reshape  the  root  locus  to 
achieve  the  desired  performance  specifications.  Generally,  three  performance' 
specifications  are  changed  by  compensation:  degree  of  stability,  transient 
response,  and  steacfy-state  error. 

Actual  condensation  is  achieved  by  addition  of  electrical  network,  and/or 
mechanical  devices  vdiich  may  contain  levers,  springs,  dashpots,  gyros,  etc. 
Biere  are  two  positions  in  the  control  system  where  condensation  is  usuedly 
performed.  In  the  feedback  loop  where  it  is  referred  to  as  feedback 
compensation  and  in  the  forward  loop  where  it  is  cedled  cascade  ccmpensation. 
In  the  forward  loop  the  compensator  is  normally  placed  in  the  low  energy  point 
so  the  power  dissipation  will  be  small. 

Common  networks  used  to  achieve  conpensation  are  lag,  lead,  and  lead-lag, 
all  of  tdiich  are  passive.  Modern  control  systems  often  use  active  networks 
which  modify  the  systems  to  ensure  the  desired  specifications  are  met 
cancelling  the  undesirable  characteristics  and  replacing  them  with  the  desired 
characteristics . 

Our  approach  to  compensation  will  be  to  look  at  our  basic  system  then  to 
see  the  effect  of  each  type  of  conpensation  on  the  shape  of  the  root  locus  and 
the  steadi-state  performance. 

13.12.1  Feedback  Compensation 

13.12.1.1  Proportional  Control  (Unity  Feedback), 
in  Figure  13.63. 


Consider  the  basic  system 


FIGUPE  13.63.  BASIC  SYSTEM 


\Aiere  for  the  first  case 


A  =  1 


G(s)  = 


y/I 


Is2  .  bs 


H(s)  =  1 


The  closed-loop  transfer  function  is 


fl3.90 
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and 

“n  '  Vij/I 

and 


ligxnre  13.64  is  the  rc»t  locus  of  this  basic  system. 


FIGURE  13.64. 

^plication  of  a  step,  R(s) 

C(s)  = 


ROOT  IDCUS  OF  BASIC  SYSTEM 
WITH  UNITY  FEEDB?.CK 


=  10/s,  results  in 


Uie  inverse  transform  of  Equation  13.91  is 


c(t) 


10  - 


-^(Ont 

e  sin  (to^t  +  (})) 
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After  a  finite  time  the  transients  die  out  and  the  re^nse  settles  down 
to  the  cormanded  value  of  10  units  with  no  steac^-state  error.  If  the  danping 

term,  b,  was  zero,  the  locus  vould  move  ri^t  to  the  imaginary  axis  and  pure 
harmonic  motion  %«Duld  result.  By  looking  at  the  error  coefficients  K^,  and 
K  ,  %je  see  that  this  ^stem  has  zero  error  for  a  stqp,  Rb/y  error  for  a  ramp 

Si 

and  has  infinite  error  for  a  parabolic  input. 

The  amount  of  viscous  danping  in  a  practical  system  is  often  limited  by 
physical  constraints.  To  oope  with  a  lightly  damped  ^^stem,  artificial 
danping  is  added.  An  investigation  of  Equation  13.90  reveals  that  the  damping 
of  the  ^tem  can  be  increased  by  increasing  the  coefficient  of  the  s  term  (e 
term)  in  the  denominator. 

13.12.1.2  Derivative  Control  (Rate  Feedback).  The  problem  of  adding  a  KC 
(Ks  term)  into  the  forward  loop  is  solved  by  derivative  control  or  rate 
feedback.  Figure  13.65  shows  the  introduction  of  the  KC  term  into  the  block 
diagram. 


T=  IC  +  bC 


FIGURE  13.65.  BASIC  SYSTEM  WITH  IKTRODUCTIC*! 

OF  KC  TEIM 

The  KC  term  can  be  achieved  by  taking  the  derivative  of  C  and  is  often 
achieved  by  means  of  a  rate  gyro  as  shown  in  Figure  13.66. 
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FIGURE  13.67.  REDUCTION  OF  FIGURE  13.66 

Ihe  control  ratio  is  written  for  Figure  13.67 

_ y _ 

Is^  +  bs  +  V  (1  +  CpS) 


C(s)  _ 
R(s) 
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y/i _ 


Cie  nat-iiral  freguency  of  the  system  is  unchanged;  however  the  danping  has 
been  increased  by  vC^/I. 

The  effect  of  derivative  control  on  the  iroot  locus  is  illustrated  in 
Figure  13.68.  The  forvard  transfer  function  becomes 

and  the  feecS^ack  transfer  function  is 

H(s)  =  (1  +  CpS) 

This  in  effect  adds  an  open-loop  zero  to  our  open-loop  transfer  function. 


GH(s)  = 


_  f  (1  ^  I 

s  (s  +  I)  ®  f ) 
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Cd(*+1/Co)  K  (8  41/Cd) 


S(8+b)  8(8  +  b/l} 

7 


WHERE  K' 


CASE  (a) 

1  <  £ 
Co  « 


CASE  (b) 


FIGURE  13.68.  RXrr  LOCUS  CF  BASIC  SYSTEM 
WITH  DERIVATIVE  OCKTRQL 


The  vctlue  of  Cp  will  determine  if  the  response  will  have  an  oscillatory 
term.  In  case  (a)  of  Figure  13.68,  1/Cj^  <  b/I  will  result  in  two  first  order 
roots  and  oscillation  is  not  possible.  However  in  case  (b)  where  l/C^^  >  b/I, 
the  locus  path  allows  for  an  oscillatory  re^nse  over  a  wide  range  of  static 
loop  sensitivity  (y/I  C^) .  It  is  (±)vious  frcm  Figure  13.68  that  high  danping 
can  be  achieved  with  reasonably  hi^  values  of  static  loop  sensitivity  (gain) . 
Before  we  can  check  the  steat^-state  error,  we  must  rearrange  the  block 
diagram  (Figure  13.69  a,  b,  c)  to  unity  feedback. 


3.12' 


m 


C(t) 

■  —  -  ■  -  -— 

- 

It*  +  bt 

1  +Cot 

(•) 


(b) 


FIGUEE  13.69.  REDOCTICN  OF  DERIVATIVE  CCKIRDL 
BLCCK  DIAGRAM  TO  UNITY  FEEDBACK 
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The  corresponding  error  coefficients  are 


K  =  «*> 

P 

K  =  ^ 

\  (b  +  uCp) 


Therefore  we  see  the  steai^— state  error  vd.ll  increase  vd.th  the 
introduction  of  derivative  control  (system  t^pe  remains  the  same) . 

We  have  discussed  two  types  of  feedback  coipensation  and  have  found  that 
the  transient  response  can  be  irtproved  by  addition  of  derivative  control,  but 
only  at  the  expense  of  the  steacfy-state  eanror.  The  other  ccrpensation 
techniques  to  be  ocaisidered  are  xased  in  the  forward  loop. 

13.12.2  Cascade  Ocaipensation 

The  first  forward  ccnpensation  technique  to  be  discussed  is  error  rate 
ccnpensation. 

13.12.2.1  Error  Rate  Compensation.  Error  rate  or  ideal  derivative 
ccnpensation  is  used  vdien  the  transient  response  of  the  systan  must  be 
ittproved.  This  is  achieved  by  reshaping  the  root  locus  by  moving  it  to  the 
left.  This  in  effect  decreases  the  system  time  constant,  t  =  therdsy 

speeding  up  the  re^xjnse.  Error  rate  coipensation  is  achieved  by  adding  the 
rate  of  change  of  the  error  signal  to  the  error  signal  (Figure  13.70) . 
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1 


FIGURE  13.70.  IDEM.  ERROR  RATE  OMPENSATOR 

The  pdiysical  effect  of  error  rate  can  be  described  as  introd\icing  anticipation 
into  the  ^stem.  The  system  reacts  not  only  to  magnitude  of  the  error,  but, 
also  to  its  probable  value  in  the  future.  If  the  error  is  changing  rapidly, 
then  E^(s)  is  large  and  the  system  responds  faster.  The  net  result  is  to 
speed  up  the  system.  Figure  13.71  shows  error  rate  added  to  the  basic  system. 


FIGURE  13.71.  BASIC  SYSTEM  WITH 
ERROR  RATE  OCNTROL 


13.130 


Now  the  forward  transfer  fimction  becomes 


P  (1  +  C  s) 
G(s)  =  — =— — — 
Is^^  +  bs 


and 

H(s)  =  1 

Therefore  the  open-loop  transfer  function  is 


We  see  that  error  rate  has  added  a  zero  to  the  function  as  well  as 
achieving  artificial  dairping  similar  to  derivative  control.  Althou^  the 
open- loop  transfer  function,  GH(s),  of  the  derivative  control  and  error  rate 
control  are  similar,  the  closed-loop  transfer  functions  are  different  due  to 
the  zero  of  the  error  rate  conpensator.  The  root  locus  for  error  rate 
is  illustrated  in  Figure  13.72.  The  zero,  in  effect,  draws  the  locus  to  the 
left;  therday  speeding  up  the  response  and  making  the  system  more  stable. 
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FIGORE  13.72.  BOOT  ICCUS  OF  BASIC  SYSTEM  WITH 
ERFDR  KATE  COMPENSATION 


Looking  at  the  error  coefficients  we  find 


K  =  ^ 
V  D 


and 

=  0 

Therefore  the  stea^^-state  performance  of  the  system  is  unchanged  frcrr*  the 
basic  system. 

In  the  real  world  an  ideal  differentiator  is  difficult  to  construct  and 
other  problens  with  differentiation  of  system  noise  arise.  A  passive  element 
known  as  a  lead  ccti^)ensator  is  used  to  approximate  ideal  error  rate  control. 
The  transfer  hmction  of  this  device  is 
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Hie  pole  tPq,  is  located  far  to  the  left  so  that  the  angle  of  the 

ocRpensator  is  nearly  all  lead  due  to  the  zero,  z  ,  vAiich  is  placed  by  trial 

o 

and  error  to  a  point  near  the  original  locus.  Hiis  normally  results  in  a 
small  increase  in  gain  and  a  large  increase  in  the  undanped  natural  frequency, 
ther^y  reducing  settling  time.  Introduction  of  this  lead  oonpensation  into 
the  basic  system  could  result  in  a  locus  as  shown  in  Figure  13.73. 


FIGURE  13.73.  LEAD  COMPENSATION  APPLIED 
TO  BASIC  SYSTEM 
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13.12.2.2  Integral  Control.  Ihe  second  cascade  carpensation  technique  is 
integral  control.  Often  the  transient  re^nse  of  a  system  is  satisfactory 
but  the  steaify-state  error  is  excessive.  integral  ccxitrol  produces  an 
actuating  signal  that  is  proportional  to  both  the  magnitude  and  the  integral 
of  the  error  signal  E(s),  Figure  13.74. 


B 


FIGUKE  13.74.  IDEAL  INTEGRAL  CCNTROL 

The  net  result  is  the  increase  in  system  type.  The  error  E^  (s)  continues 
to  increase  as  long  as  an  error,  E(s) ,  is  present  and  eventually  becanes  large 
enou^  to  produce  an  output  signal  equal  to  the  iipfut.  Uie  error,  E(s),  is 
then  zero.  Since  we  are  not  interested  in  changing  the  time  response  of  the 
^tem,  the  positioning  of  the  zero  becanes  very  important.  The  pole  at  the 
origin  has  the  effect  of  moving  the  locus  to  the  right  and  slewing  down  the 


effect  of  the  pole  on  the  locus. 

Figure  13.75  shews  integral  control  added  to  out  basic  system. 


FIGURE  13.75.  BASIC  SYSTEM  WITH  INTEGRAL  CONTROL 


New  tlie  forward  transfer  function  beccnces 


G(s)  = 


VI  (s  +  C^) 
s(Is^  +  bs) 


i  (s  +  c. ) 


H(s)  =  1 


Ttie  open-loop  transfer  function  is 


GH(s)  = 


i  (s  +  C . ) 


(■•I) 


(13.92) 


Figure  13.76  is  the  root  locus  of  Equation  13.92. 
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LIES  IN  RHP  AND  SYSTEM  IS 
UNSTABLE  FOR  ALL  K) 

FIGURE  13.76.  ROOT  LOCUS  OF  BASIC  SXSTFM  WITH 
INTBaRAL  OOmOL 

Die  pole  at  the  origin  increases  the  system  type  as  well  as 

the  Cf^aiing  of  K  along  the  loci.  The  change  in  the  transient  re^nse  for  the 
ideal  integral  control  is  minimized  by  placing  the  zero  very  close  to  origin. 
In  a  real  world  system  there  is  a  limit  to  achieving  this  prcKimity. 

The  integral  control  oonpensator  is  achieved  electrically  by  a  lag 
network  or  mechanically  by  lase  of  an  integrating  gyroscope. 

The  lag  network  transfer  function  is  of  the  form 


vAiere  a  >  1  and  is  usually  about  10;  ther^ore.  Equation  13.93  is 
approKunately 

G  (s)  s  A  (s  +  1/t) 
c  ot  s 

In  design,  if  the  pole  and  zero  are  placed  very  close  (at  the  origin), 
the  net  angle  contribution  at  the  dondnate  poles  can  be  kept  to  less  than  5° 
and  the  locus  is  displaced  only  sli^tly.  Ihe  resulting  increased  gain  of  the 
^stem  and  increased  ^stem  type  all  result  in  decreased  steac^-state  error. 

The  error  coefficient  for  the  ccrapensated  system  became 

Kp  =  « 


\  = 


and  the  system  will  now  handle  a  parabolic  input  with  a  e(t)  =  R  b/jiC. . 

SS  1 

A  suiraiary  of  passive  conpensations  contained  in  Table  13.5. 
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PHYSICAL  SYSTEM 


13.13  SUMMARY 


In  this  chapter  an  attempt  has  been  made  to  present  the  fundanentals  of 
control  system  analysis,  implications  of  the  theory  were  held  to  a  minimum  so 
that  full  attention  could  be  devoted  to  learning  the  tools  and  techniques  used 
in  this  of  analysis. 

Once  the  analysis  techniques  have  been  mastered,  the  more  interesting 
appropriate  areas  of  fli^t  control  and  handling  qualities  may  be  addressed. 
A  knowledge  of  root  locus  theory  and  frequency  response  is  essential  in 
understanding  the  aiplications  of  feedback  analysis  to  fli^t  vehicle  systems. 

Despite  the  introduction  of  modem  control  theory  (state  variables)  aid 
digital  flight  control,  an  understanding  of  these  systems  is  still  based  in 
large  measure  on  knowledge  of  classical  feedback  control  systems. 
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